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FOREWORD 


This  study  of  water  entry  impact  shock  on  flat  faced  cylindrical  missiles 
was  sponsored  by  Dr.  Thomas  Peirce  of  NAVSEA,  Code  63R31,  whose  continued 
support  of  water-entry  research  has  made  this  extension  of  prior  work  possible. 

This  work  may  be  applied  to  the  determination  of  the  shock  component  of 
water  entry  impact  force,  as  a  function  of  time,  on  air  dropped  torpedos  and 
mines  with  stabilizing  flat  noses.  The  program  IMPACT  lessens  the  need  for  the 
much  more  expensive  water  entry  impact  testing  on  prototypes  or  models. 


Approved  by: 


WILLIAM  C.  JONES 
Surface  ASW  Division 


."'U  . 

a,  Acting  Head 


i  Aoconalon  For  ! 

I'ntis 

CIRA&I 

. XT 

□ 

DTIC  TAB 

Unannounced 

1:1 

'  Just  If  lent,  ion .  .J 

By.,  .. 
Din  •• 

....  .  .. 

. — -- 

1  s  m/ 

Av  .1 

. .  i . >  i  ■ 

A  .':i  i  l  .‘i 1 1 •  J /  u  V  | 

lal 

PHSCiDlNa  PAOI  BLANK-, NOT  AIMED 


1 


NSWC  TR  82-438 


PREFACE 


The  junior  author  would  like  to  express  his  appreciation  to  Mr.  Charles  W, 
Smith  and  Dr.  John  L.  Baldwin  of  NSWC 1 s  Hydroballistic  Facility,  Code  U23 ,  for 
their  direction  and  support,  and  in  making  this  study  available  as  a  rotational 
assignment  project. 


it 


NSWC  TR  82-438 


CONTENTS 

■Lage 

INTRODUCTION . . .  1 

PROBLEM  DESCRIPTION  .  1 

GEOMETRY  OF  IMPACT .  1 

PHYSICS  OF  IMPACT  SHOCKS .  3 

BASIC  ASSUMPTIONS  FOR  SOLUTION .  6 

DERIVATION  OF  IMPACT  SHOCK  DYNAMICS,  .  8 

MISSILE  FACE  DEFORMATION  RESPONSE  .  8 

WATER  PARTICLE  VELOCITY  .  14 

RANKINE-HUGONIOT  SHOCK  PRESSURE  .  20 

EXPANSION  WAVE  SPEED  ACROSS  DEFORMING  MISSILE  FACE . 25 

TIME  DURATION  OF  IMPACT  SHOCK . 28 

SHOCK  AREA  ON  MISSILE  FACE.  .  . . 29 

IMPACT  SHOCK  FORCE  AND  IMPULSE.  '.  .  .  . . 39 

IMPACT  SHOCK  DYNAMICS  COMPUTATION  PROGRAM:  'IMPACT' . 40 

DESCRIPTION  OF  PROGRAM  STRUCTURE . 40 

USE  OF  THE  PROGRAM . 43 

EXAMPLE  PROBLEM  .....  .  44 

RESULTS  OF  THE  EXAMPLE . 45 

POSSIBLE  MODIFICATIONS  AND  IMPROVEMENTS  TO  THE  PRESENT  FORMULATION  .  48 

NOMENCLATURE . 49 

Appendix  Page 


A  'IMPACT'  PROGRAM  LISTING  .  A-l 

B  LISTING  OF  EXAMPLE  OUTPUT .  B-l 


ill 


NSWC  TR  82-438 


ILLUSTRATIONS 


;ure  Page 

1  GEOMETRY  OF  A  SHOCK  PRODUCING  OBLIQUE  ENTRY  IMPACT  .  A 

2  MISSILE  FACE  DEFORMATION  UNDER  SHOCK  PRESSURE  AREA  .  7 

3  WAVE  TRANSMISSION  ALONG  A  UNIFORM  BAR  WITH  NO 

LATERAL  OR  TRANSVERSE  DISPLACEMENTS . 10 

4  MISSILE  FACE  DEFORMATION  AND  THE  WATER  SHOCK  WAVE  FRONT . 15 

5  WATER  PARTICLE  VELOCITY  IN  THE  SHOCKED  REGION . 16 

6  OBLIQUE  SHOCK  WAVE  GEOMETRY . 19 

7  SHOCK  WAVE  SPEED  IN  WATER  AS  A  FUNCTION  OF  SHOCK  PRESSURE . 22 

8  SECANT  ITERATION  TECHNIQUE  .  24 

9  GEOMETRY  OF  EXPANSION  WAVE  MOVING  ACROSS 

DEFORMING  MISSILE  FACE  .  26 

10  WATER  ENTRY  AND  IMPACT  SHOCK  DURATION  TIMES . 30 

11  SHOCK  AREA  ENVELOPE  DURING  WETTING  ENTRY  .  32 

12  SHOCK  AREA  ENVELOPE  AFTER  FULL  WATER  ENTRY  .  36 

13  SIMPSON'S  RULE  INTEGRATION  FOR  A  PARAMETRIC  CURVE . 3? 

14  IMPACT  SHOCK  FORCE  AS  A  FUNCTION  OF  WATER  ENTRY  TIME  .  45 

15  MAXIMUM  IMPACT  SHOCK  FORCE  AS  A  FUNCTION  OF  ENTRY  ANGLE . 46 

16  TOTAL  IMPULSE  AS  A  FUNCTION  OF  ENTRY  ANGLE  .  47 


iv 


NSWC  TR  82-438 


INTRODUCTION 


This  study  is  an  extension  of  a  paper  by  Dallas  D.  Laumbach^  of  water 
entry  impact  shock  on  flat  faced  cylindrical  missiles  with  their  axis  normal  to 
the  water's  surface.  This  extension  includes  shock  producing  oblique  entry 
impacts  as  well  as  normal  impacts.  The  results  of  this  work  are  such  that  for 
the  normal  entry  case  they  reduce  to  those  of  I.aumbach's. 


PROBLEM  DESCRIPTION 

The  problem  treated  by  this  report  is  the  determination  of  the  response 
dynamics  of  a  flat  faced  cylindrical  missile  striking  a  smooth  water  surface  in 
a  shock  producing  impact.  The  missile  enters  the  water  with  its  axis  parallel 
to  its  velocity  vector  and  either  normal  or  oblique  to  the  water's  surface.  The 
dynamic  response  of  a  missile  to  water  entry  impact  shock  can  be  described  by 
the  impact  shock  force  acting  on  the  missile's  face  as  a  functipn  of  water  entry 
time  and  its  total  impulse.  To  determine  these  quantities  it  is  necessary  to 
find  the  shock  pressure,  resulting  from  impact,  between  the  water  and  the 
missile's  face.  This  requires  consideration  of  shock  wave  propagation,  compres¬ 
sibility  of  water  in  shock,  deformation  of  the  missile's  face  under  shock 
pressure,  water  particle  speeds,  and  the  geometry  of  the  relieving  expansion 
wave.  Only  the  force  due  to  shock  pressure  is  found;  gravity,  buoyancy,  and 
drag  forces  were  not  treated.  The  general  problem  for  oblique  entry  is 
intrinsically  three  dimensional  with  time  varying  boundary  conditions,  due  to 
gross  penetration,  coupling  the  compressing  water  and  deforming  missile. 


GEOMETRY  OF  IMPACT 

Described  below  is  both  the  geometry  of  impact  defining  the  problem  to  be 
solved  and  a  nomenclature  of  phrases  which  will  be  used  in  reference  to  this 
geometry  throughout  the  report. 

1)  The  missile  is  a  circular  cylinder  with  radius  r  and  of  semi-infinite 
length . 


^Laumbach,  D.D.,  Impact  Shock  on  a  Blunt-Nosed  Missile  Striking  the  Water  with 
its  Axis  Normal  to  the  Water  Surface,  ^andia  Corp.  SC-TM~64"987 ,  Sept  1964. 
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2)  The  end  of  this  cylinder  which  strikes  the  water  is  initially  a  plane 
surface  normal  to  the  missile's  axis.  This  is  called  the  missile's  face  and  its 
undeformed  configuration  is  referred  to  as  the  missile's  rigid  face. 

3)  The  undisturbed  water  is  a  semi-infinite  half  space  bounded  by  a  plane 
called  the  water's  surface.  This  undisturbed  water  has  an  initial  density  of 

p  i  at  a  pressure  of  pi,  1  atmosphere,  and  a  temperature  of  20°C. 

4)  The  missile  enters  the  water  with  a  constant,  and  prescribed,  velocity 
vector  which  is  parallel  to  the  missile's  axis. 

5)  This  velocity  vector  is  prescribed  by  a  speed  V'  and  an  orientation 
a,  the  angle  between  the  missile's  axis  and  a  normal  to  the  water's  surface. 

6)  Water  entry  impact  time,  t,  is  the  time  since  the  missile  first 
contacts  the  water's  surface. 

7)  For  oblique  water  entry,  a>0,  the  point  which  first  touches  the 
water's  surface  is  called  the  initial  contact  point  and  its  associated  generator 
of  the  missile's  cylinder  is  called  the  initial  entry  edge. 

8)  Also  for  oblique  entry,  the  contact  between  the  missile's  face  and  the 
undisturbed  water  surface  is,  for  t>0,  a  chord  on  the  missile's  face  called 

the  wetting  line. 

9)  Note  that  for  normal  water  entry,  a^O,  the  full  missile  face  is  in 
initial  contact  with  the  water's  surface  at  t"0  and  there  is  no  wetting  line. 

10)  The  wetting  line  moves  across  the  undisturbed  water  surface  with  the 
speed  V'/sina,  called  the  surface  disturbance  speed. 

11)  This  surface  disturbance,  caused  by  the  moving,  wetting  line,  results 
in  a  discontinuous  disturbance  of  the  water  particles,  imparting  a  velocity,  and 
raising  the  water's  density  and  pressure  to  p f  and  pf  respectively.  This 
disturbance,  a  shock,  propagates  radially  as  rays  from  the  instantaneous  wetting 
line  position  with  a  speed  a,  the  shock  wave  speed  in  water. 

12)  These  water  shock  rays  may  form  the  loci  of  a  plane  surface  which 
subtends  the  wetting  line  and  moves  through  the  water  with  speed  a.  This 
surface  is  the  shock  wave  front. 

13)  The  wetting  line  moves  across  the  face  of  the  missile  with  a  speed 
V“V'/tan  o,  called  the  wetting  line  speed. 

14)  Water  in  a  shocked  state  recovers  its  initial  properties  in  an 
expansion  which  propagates  through  the  shocked  water  as  a  wave,  moving  with 
speed  c',  forming  the  expansion  wave  front. 

15)  The  contour  of  this  expansion  wave  front  on  the  missile's  face  is 
called  the  expansion  wave  profile  and  moves  over  the  missile's  face  with  a 
projected  speed  c  in  the  rigid  face. 
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16)  The  region  of  the  missile's  face  between  the  wetting  line  and  Che 
expansion  wave  profile  encompasses  water  in  shock,  deforming  the  missile's  face 
with  shock  pressure  pf.  This  region  is  called  the  shock  pressure  region  and 
its  area  the  shock  pressure  area  S. 

17)  The  missile's  entry  speed  and  angle  are  such  that  its  impact  is  shock 
producing.  (This  requires  that  the  surface  disturbance  speed  be  greater  than 
the  water's  shock  wave  speed,  establishing  the  shock  wave  front  at  tho  wetting 
line,  and  further,  that  the  wetting  line  speed  is  greater  than  the  expansion 
wave  profile  speed,  resulting  in  a  shock  pressure  region.) 

18)  The  plane  which  contains  the  missile's  axis  and  is  normal  to  the 
water's  surface  is  called  the  plane  of  symmetry.  This  plane  contains  the 
initial  contact  point,  bisects  the  wetting  line,  and  is  a  symmetry  plane  to  the 
expansion  wave  profile  and  shock  pressure  area. 

19)  The  loading  by  the  water  shock  pressure  on  the  missile's  face 
initiates  a  deformation  of  the  face  at  the  wetting  line  which  propagates  strain 
into  the  missile  as  a  shock  wave  as  the  shock  pressure  region  moves  across  the 
missile's  face.  This  missile  shock  wave  has  an  axial  speed  of  aro. 

20)  The  water  expansion  wave,  restoring  the  water  to  its  initial  pressure, 
unloads  the  missile's  deforming  face  of  its  shock  pressure.  This  unloading 
propagates  its  effects  as  a  wave  restoring  the  missile's  material  to  its  initial 
undisturbed  state. 

The  geometry  of  missile  impact  for  oblique  entry  is  illustrated  by  Figure 
1,  which  depicts  the  plane  of  symmetry  and  the  axial  projection  of  the  missile's 
face. 


PHYSICS  OF  IMPACT  SHOCKS 

The  physics  of  shocks  in  water,  as  distinct  from  gasses,  is  treated  in 
depth  by  Robert  H.  Cole,  with  regard  to  underwater  explosions,  in  Chapters  2  and 
4  of  his  work. ^  The  impact  shock  response  of  solids  is  summarized  in  the 
first  and  second  chapters  of  the  reference  work  by  J.A.  Zukas,  et  al3  which 
deals  with  the  recent  techniques  of  analysis,  including  numerical,  and  experi¬ 
mental  impact  dynamics.  What  follows  is  e  brief  sketch  of  the  response  physics 
of  the  water  and  missile  to  impact  shock  as  it  applies  to  the  problem  of  this 
study. 


^Cole,  R.H.,  Underwater  Explosions  (Princeton:  Princeton  University  Press, 
1948). 

^Zukas,  J.A.,  Nicholas,  T.,  Swift,  H.F.,  Greszczuk,  L.B.,  and  Curran,  D.R. , 
Impact  Dynamics  (New  York:  John  Wiley  &  Sons,  Inc.,  1982). 
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FIGURE  1 .  GEOMETRY  OF  A  SHOCK  PRODUCING  OBLIQUE  ENTRY  IMPACT 
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Water  Impact  Response 

Water  entry  by  the  missile  imparts  a  disturbance  speed  to  the  water 
particles  in  contact  with  the  missile's  wetting  face,  initiating  at  the  wetting 
line.  These  particles  propagate  their  disturbance  to  the  rest  of  the  water  as  a 
compressive 

shock  wave,  forming  the  discontinuity  between  disturbed  and  undisturbed  regions 
of  water.  The  rays  of  this  shock  wave  are  generated  by  the  movin.v  wetting  line 
and  propagate  at  the  shock  wave  speed,  a,  from  the  position  the  wetting  line  had 
when 

the  ray  was  first  generated.  With  the  wetting  line  moving  across  the  water's 
surface  faster  than  the  shock  wAve  speed,  the  rays  it  generates  form  the  loci  of 
a  plane,  the  shock  wave  front,  which  subtends  the  wetting  line  At  an  angle,  the 
Mach  angle,  of  less  than  90°.  With  the  missile  then  penetrating  the  shocked 
water,  the  water  is  compressed,  resulting  in  an  increase  in  density  pf  and 
pressure  pf,  the  shock  pressure . 

Along  the  wetting  circumference  of  the  missile  the  water  is  free  to  expand 
and  recover  its  initial  state.  This  expansion  propagates  as  a  wave,  the 
expansion  wave,  whose  rays  are  generated  at  t^he  instantaneous  positions  of  the 
cusps  of  the  wetting  line  on  the  missile's  circumference.  These  rays  form  the 
loci  of  an  expansion  wave  surface  chasing  the  shock  wave  front.  The  water  being 
compressed  with  the  passage  of  the  shock  wave  front  is  then  later  restored  to 
its  initial  state  with  the  passage  of  this  expansion  wave.  Since  water  is 
unable  to  sustain  tensile  stresses  this  expansion  wave  is  weak  compared  with  the 
shock  wave  and  moves  through  the  water  at  the  sonic  speed  c',  which  is  L®.bs  than 
the  shock  wave  speed.  This  trailing  of  the  expansion  wave  results  in  a  shock 
pressure  region  contained  between  the  shock  wave  front  and  the  expansion  wave 
surface.  On  the  face  of  the  missile  the  shock  region  is  contained  between  the 
wetting  line  and  the  expansion  wave  profile  on  the  missile's  deforming  face. 

For  the  special  case  of  normal  water  entry  the  face  of  the  missile  is  in 
full  contact  with  the  water's  surface  at  the  instant  of  impact,  generating  a 
shock  wave  front  parallel  to  the  missile's  face  with  a  Mach  angle  of  zero. 

Also,  the  expansion  wave  is  generated  simultaneously  on  the  full  circumference 
of  the  missile's  face  resulting  in  a  circular  expansion  wave  profile  concentric 
with  the  circumference,  and  whose  radius  diminishes  with  speed  c. 


Missile  Impact  Response 

For  oblique  water  entry,  shock  pressure  acting  on  the  missile's  face  first 
loads  the  face  at  the  wetting  line  as  it  moves  across  the  face  with  speed  V. 

This  loading  results  in  compressive  axial  deformation  of  the  face  and  the 
propagation  of  a  compressive  shock  initiating  at  the  wetting  line.  This  shock 
wave  moves  up  the  length  of  the  missile  as  axial  rays  generated  by  the  instanta¬ 
neous  position  of  the  wetting  line  and  travel  with  the  wave  speed  a,,,,  forming 
an  oblique  wave  front.  The  face  continues  to  deform  obliquely  from  the  wetting 
line  in  the  moving  shock  pressure  region  with  an  assumed  constant  deformation 
angle,  w,  away  from  the  rigid  face  position.  This  region  of  the  face  is  then 
unloaded  with  the  passage  of  the  water's  expansion  wave  profile,  moving  across 
the  face  at  speed  c,  at  which  time  compressive  deformation  ceases.  This 
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unloading  also  propagates  its  effects  as  an  expansion  wave  which,  since  metals 
sustain  full  reversed  tensile  stresses,  also  travels  at  the  shock  wave  speed 
am.  This  expansion  wave  in  the  missile,  being  generated  by  the  moving  water's 
expansion  wave  profile,  creates  a  second  oblique  wave  front  in  the  missile 
restoring  the  missile's  material  to  its  initial  unstrained  state.  This  process 
of  missile  face  deformation,  shock  wave  propagation,  and  expansion  wave 
unloading  is  illustrated,  for  the  stages  of  water  entry  time  Te  and  impact 
shock  duration  time  T^,  by  Figure  2,  which  shows  the  missile  in  the  plane  of 
symmetry . 

For  the  case  of  normal  water  entry  impact,  with  the  full  face  exposed  to 
shock  pressure  at  the  instant  of  impact,  the  shock  wave  front  and  deformation 
surface  become  planes  normal  to  the  missile's  axis,  the  deformation  nngle  is  is 
zero.  Also  with  the  water's  expansion  wave  moving  in  from  the  full  circumference 
simultaneously,  unloading  the  shock  pressure  region  at  speed  c,  the  deformation 
profile  in  all  axial  planes  becomes  a  truncated  cone  whose  top  progresses  upward 
to  form  a  full  cone  at  shock  duration  time  T^. 


BASIC  ASSUMPTIONS  FOR  SOLUTION 

Because  a  three  dimensional  elasticity/non-ideal  compressible  fluid 
dynamics  model  with  variable  boundary  conditions  is  so  complex,  even  for  a 
finite  difference  approach  on  the  resulting  system  of  coupled  nonlinear  partial 
differential  equations,  a  simplified,  essentially  one  dimensional,  "engineering" 
type  of  solution  is  sought.  Such  a  solution  requires  a  broad  base  of  simplifying 
assumptions  which  reduce  the  problem  to  a  soluble  model  without  sacrificing  the 
essential  geometry  and  physics  of  the  real  phenomena.  The  basic  assumptions 
employed  for  this  solution  are  enumerated  below: 

1)  The  missile's  water  entry  speed  V'  is  constant  during  impact  shock. 

2)  The  missile's  water  entry  angle  a  is  also  constant  during  impact 
shock . 

3)  Water  shock  is  an  adiabatic,  inviscid,  process  with  shock  pressure  and 
wave  speed  empirically  related. 

4)  The  impact  shock  pressure  p£  is  constant  in  time  and  uniform  over  the 
shock  pressure  region  on  the  missile's  face. 

5)  The  missile's  face  under  shock  pressure  deforms  with  a  constant  axial 
speed  w,  initiated  at  the  wetting  line  and  ceasing  at  the  expansion  wave 
profile,  resulting  in  a  planar  deformation  surface  spanning  the  shock  pressure 
region.  This  deformation  is  described  by  an  axial  displacement  w,  which  is  the 
‘projection  of  the  deformation  plane  to  the  initial  entry  edge  and  is  measured 
along  this  edge  from  the  undeformed  face  position  as  shown  in  Figure  1. 

6)  This  deformation  is  small  euouqh  so  as  to  result  in  a  linear  strain- 
dioplacement  relation. 
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7)  The  missile  is  a  homogeneous,  isotropic,  linearly  elastic  material  o£ 
sufficient  length  so  as  to  preclude  shock  wave  interaction  by  reflections  during 
impact  shock.  The  missile  is  thus  characterized  by  its  initial  or  undisturbed 
properties  of  density  pm,  Young's  modulus  E,  ai.  i  Poission's  ratio  v. 

8)  The  axial  cylinder  of  the  shock  pressure  region  in  the  missile  is 
constrained  to  have  no  lateral  or  transverse  strains. 

9)  Both  compression  and  expansion  shock  waves  in  the  missile  are  one 
dimensional,  moving  axially  with  a  speed  am.  The  compression  wave  is 
initiated  at  the  wetting  line  and  the  expansion  wave  is  initiated  at  the  water 
expansion  wave  profile. 

10)  The  water  expansion  wave  is  sonic,  moving  with  the  speed  of  sound  c1 
in  water  and  restoring  the  water  to  its  initial  pressure  p£. 

With  the  problem  defined  and  the  assumptions  established  a  derivation  of 
impact  shock  dynamics  is  now  carried  out,  leading  to  the  determination  of  shock 
force  and  impulse. 


DERIVATION  OF  IMPACT  SHOCK  DYNAMICS 


The  equations  which  are  the  basis  for  determining  the  impact  shock  force 
and  total  shock  impulse  arc  now  derived.  This  derivation  is  based  upon  the 
geometry,  physics,  and  simplifying  assumptions  which  were  established  in  the 
Introduction.  These  allow  an  elementary  or  "engineering"  solution  of  an 
otherwise  complex  problem  in  mechanics  which  it  hopefully  approximates. 


MISSILE  FACE  DEFORMATION  RESPONSE 

As  the  missile  enters  the  water,  impact  shock  pressure  produces  a 
compressive 

deformation  over  the  shock  area  of  the  wetted  face.  This  axial  displacement  has 
been  assumed  to  progress  elastically  with  constant  particle  speed  w(z,t),  at  the 
time  of  wetting,  resulting  in  the  linear  deformation  profiles  in  the  plane  of 
symmetry  as  depicted  by  Figure  2.  Thus  the  displacement  w  of  the  deforming  face 
extended  to  the  initial  entry  edge,  as  shown  in  Figure  1,  is 


w(z”0, t )  *  w(z*0,t)t . 


The  particle  speed  may  be  expressed  as 


w(t,z)  -±”  I* 
3z  dt 
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where 


8w 
3  z 


Lim 
Az  + 


[<w 

0 


-  Ai i  ’Az)  -  w]/Az 
Az 


is  the  change  in  the  length  of  Az  per  unit  length  caused  by  the  displacement 
field  w(z,t)  acting  along  Az,  i.e.  the  axial  strain 


and 


AlS.  ■  Lim  AS. 

dt  At 

At  0 


is  the  distance  Az  over  which  the  particles  have  been  displaced  per  unit  time 
At,  i.e.  the  speed  with  which  the  particles  have  been  disturbed  -  the  shock 
wave  speed  in  the  missile 


dz 

dt 


Therefore  the  displacement  at  the  initial  contact  point  becomes 

w  -  CjsBnt,  (1) 

where  c2  ■  cz(z“0,t)  is  the  strain  at  the  initial  contact  with  t  being 
the  time  since  initial  water  contact  by  the  missile. 


Compressive  Elastic  Wave  Speed 

The  missile  has  been  assumed  to  respond  elastically  to  impact. 

Furthermore,  since  the  shock  pressure  area  on  the  missile  face  is,  except  for 
normal  impact  at  t"0,  contained  within  the  radius  of  its  sides,  and  thus 
constrained  by  them,  the  elastic  wave  speed  is  derived  as  a  one  dimensional  wave 
transmission  along  a  uniform  bar  constrained  to  have  no  lateral  or  transverse 
displacements.  This  derivation  is  given  by  W.  Johnson*  and  is  developed  below 
with  the  general  situation  illustrated  by  Figure  3. 


^Johnson,  W. ,  Impact  Strength  of  Materials  (London:  Edward  Arnold  Publishers, 
Ltd.,  1972),  pp.  13-14. 
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- ► 

w(i,t) 


FIGURE  a  WAVE  TRANSMISSION  ALONG  A  UNIFORM  EAR  WITH  NO  LATERAL  OR  TRANSVERSE  DISPLACEMENTS 
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The  Hooke's  Law  stress-strain  relation  for  linear  elastic  materials 
for  example  Timoshenko  and  Goodier,^ 

ex  ■  I  (ox  -  v(oy  +  oz)] 

Ey  -  I  [oy  -  V(cx  +  02)j 

Cjj  “  I  loz  -  v(ox  +  Oy)]  , 

which  with  the  zero  lateral  and  transverse  strain  assumption 
cx  -  ey  -  0 
becomes 

0  “  ox  -  v(0y  +  o2) 

0  -  Oy  -  v(cx  +  o2) 

c2  -  I  [o2  -  v(ox  +  Oy)]  . 

E 

These  require  that 

o  y  ■  ov  ■  -H—  c,  and 
x  y  1-v  z 


^Timoshenko ,  S.P.  and  Goodier,  J.N.,  Theory  of  Elasticity.  3rd.  ed.  (New 
McGraw-Hill  Book  Co.,  Inc.,  1970),  p.  8. 


is,  see 


(2) 


York: 
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With  the  strain-displacement  relation 


a  stress-displacement  relation  is  found  by  substitution  for  cz  as 

3w  „  (l+v)(l-2v) 

32  E  ( 1 — v ) 


or 


30  z  m  E(  1-v )  3^w 

3z  (l+v)<!l-2v)  8z2 


A  balance  of  forces  with  momentum  acting  on 
shock,  as  shown  in  Figure  3,  yields  the  equation 


an  element  of  the  material 
of  motion 


in 


o2S 


3o 

+  _ £  «z)S 

3z 


pra  (S6z) 


3tJ 


or 


3o 
_ s 

3z 


-  o 


m 


afw 

3t2 


Substitution  for  stress  oz  from  the  previously  determined  stress-displacement 
relation  gives  the  wave  equation 


3 ^w  .  Hji.  (1-v)  3 2w 

3t2  Pm  Tr+v)(l-2v)  3z2 


■  ■ 
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Which  implies  an  elastic  wave  speed  in  the  missile  of  the  form 


■v 


E 

pm 


( 1— v ) 

<l+v)"(l-2v) 


(3) 


since 


32w 


3f 


32w 

3z2 


has  the  general  solution 

w(z,t)  -  f(z-amt)  +  g(z+arat)  , 

which  are  waves  with  profile  forms  f(z)  and  g(z)  translating  along  the  z-axis  to 
the  right  and  left  respectively  with  wave  speed  am. 

Axial  Strain  Under  Shock  Pressure 

On  the  shock  pressurized  face  of  the  missile  the  axial  stress'  is  the  shock 
pressure 

o2(0,t)  -  -pf  . 

Substituting  this  into  Equation  (2)  yields  the  strain  at  the  initial  contact  point 


«-  (n  .  Pf  .  ( 1+v ) ( l-2v  ) 

ezto,t;  r  — - 


Then  introducing  the  wave  speed  from  Equation  (3)  gives  this  strain-shock  pressure 
relation  the  form 


c 


z 


-pf 

■  i  -  —  e 

Pm  am 


(4) 
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Deformation  Angle  from  a  Rigid  Face 

With  the  initial  water  entry  contact  point  displacement  found  the  angle  at 
which  the  face  deforms  from  a  rigid  one  in  the  plane  of  symmetry  may  now  be 
determined.  The  geometry  of  this  deformation  is  illustrated  by  Figure  4.  As  can 
be  seen  from  this  figure  the  deformation  angle  is 


cat 
tan  u)  " — z  — 
V1  t 

tan  a 


or 


.  /c  a  tan  a  \ 

■  t*"-1  — j  • 


WATER  PARTICLE  VELOCITY 

The  development  of  water  shock  dynamics  requires  a  knowledge  of  the  water 
particle  velocity  in  the  shock  region,  particularly  that  component  normal  to  the 
shock  wave  front.  The  first  step  is  to  determine  the  resultant  water  particle 
velocity  in  thia  region  relative  to  the  water  at  rest,  or  equivalently,  to  the 
initial  water  entry  contact  point. 


Resultant  Water  Particle  Velocity 

As  the  missile  enters  the  water  its  face  imparts  motion  to  the  water 
particles  in  the  shocked  region  which  move  with  the  speed  of  the  deforming  face 
and  in  the  direction  normal  to  it.  This  deforming  face  speed  relative  to  the 
fixed  entry  point  is  the  entry  speed  of  the  missile  less  the  deformation  speed  of 
the  missile's  face  relative  to  the  missile,  or  w  of  Equation  (1).  The  geometry  of 
these  velocity  vectors  acting  on  a  water  particle  in  the  shocked  region  is 
depicted  by  Figure  5.  The  component  of  this  difference  in  the  direction  normal  to 
the  shock  deformed  face  plane  is  the  resultant  water  particle  velocity  in  the 
shock  preasure  region, 


u  ■  (V'  -  e  zam)  cos  u> 


(6) 


Mach  Angle  of  Water  Shock  Wave  Front 

With  the  wetting  line  moving  across  the  surface  of  the  water  at  a  speed  V, 
the  surface  disturbance  speed,  greater  than  the  shock  wave  speed  in  water  a  shock 
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FIGURE  4.  MISSILE  FACE  DEFORMATION  AND  THE  WATER 
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wave  front  is  established  from  the  wetting  line.  The  angle  at  which  this  wave 
front  subtends  the  undisturbed  surface  is  the  Mach  angle.  The  geometry  is 
illustrated  in  Figure  4,  from  which 


sin  13  “  -5^ — 
V 1 1 
sin  a 


or 


(7) 


For  a  Mach  angle,  or  shock  wave  front,  to  subtend  the  wetting  line  it  is 
necessary  that 


a  sin  a  ✓ 

~T' 


1. 


Thus  the  following  cases  are  possible: 


1) 


sin  a 


<  a  impact  is  nonshock  producing;  or 


2)  ~1 —  >  a  then  impact  shock  occurs  and  either: 


sm  a 

a)  V  <  a  for  which  R  >  a, 

b)  V'  ■  a  for  which  13  ■  a,  or 

c)  V'  >  a  for  which  13  <  a 


with  the  boundary  requiring  that  13  _>  a-w,  where  13  ■  a-w  is  equivalent  to 
normal  entry  impact. 


Critical  Entry  Angle 

Equation  (7)  provides  one  necessary  condition  for  water  impact  to  be  shock 
producing,  namely,  for  a  chock  wave  front  to  be  established  it  is  necessary  that 


V1 

sin  a 


> 


a, 


(8) 
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or  the  speed  of  the  disturbance  across  the  surface  of  Che  water  must  be  greater 
than  the  shock  wave  speed  in  water.  The  other  necessary  condition  for  water 
impact  to  be  shock  producing  is  that  the  wetting  line  speed  across  the  face  of  the 
missile  must  be  greater  than  the  expansion  wave  speed  across  the  projected  rigid 
face,  or 


V 


_y'_  > 

tan  a 


c . 


(9) 


This  is  necessary  for  a  nonzero  shock  pressure  area,  as  can  be  seen  in  Figure  L. 

The  critical  entry  angle  a*  is  the  entry  angle  of  the  missile  beyond  which 
the  impact  is  no  longer  shock  producing.  From  the  two  necessary  conditiona  on 
shock  impact,  Equations  (8)  and  (9),  the  critical  entry  angle  will  be  the  smaller 
of 


si"-‘  ft) 

(10) 


\  *) 


Component  of  Water  Particle  Velocity  Normal,  to  Shock  Wave  Front 

With  the  missile  face  deformation  angle,  resultant  water  particle  velocity, 
and  water  shock  wave  front  Mach  angle  found,  the  component  of  water  particle 
velocity  normal  to  the  shock  wave  front,  q,  may  now  be  determined.  The  geometry 
of  these  water  particle  velocity  vectors  is  illustrated  by  Figure  6.  With  the 
reoultant  water  particle  velocity  u  perpendicular  to  the  missile's  deforming 
missile  face  and  the  normal  component  q  perpendicular  to  the  shock  wave  front  the 
angle  between  them,  y ,  is  the  angle  between  the  deforming  missile  face  and  the 
shock  wave  front,  or 


a* 


min 


y  ■  13  -  (a-w) 


(11) 


Thus  the  normal  component  of  water  particle  velocity  is 


q  *  u  cos  y 


(12) 
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Note  that  for  the  special  case  of  normal  water  entry  impact,  a“0,  Equations 
(5),  (6),  (7),  (11),  and  (12)  reduce  to: 

to  ■  0 , 

u  ■  v’  "  ezam> 

13  “  0, 
y  m  0,  and 
q  -  u; 

which  is  consistent  with  the  results  of  D.  Laumbach's  paper. 


RANKINE-HUGONIOT  SHOCK  PRESSURE 

A  disturbance  of  finite  strength,  a  shock,  will  propagate  through  water  as  a 
compressive  wave  with  a  speed  greater  than  the  local  speed  of  sound.  These  shocks 
move  too  fast  for  gross  temperature  changes  to  occur,  i.e.  adiabatically .  Thus 
the  relation  between  shock  wave  speed,  pressure,  and  particle  velocity  may  be 
derived  by  use  of  only  the  balance  of  momentum  and  the  continuity  equations.  This 
is  described  by  C.-S.  Yih**  and  the  derivation  of  the  shock  pressure,  or 
Rankine-Hugoniot  relation,  fallows  below. 


Derivation  of  the  Rankine-Hugoniot  Relation 

'Consider  flow  through  the  shock  wave  front  shown  in  Figure  6.  With  the  wave 
sweeping  through  a  mass  of  water  of  unit  cross  sectional  area  in  unit  time  we  have 
from  the  - 

■ 

Continuity  Equation: 


p ia  ■  p  f  (a-q)  and 


Momentum  Balance: 


Pf  "  Pi  “  P ia 


Pf 


(a-q)2 


k  Yih,  C.-S.,  Fluid  Mechanics/A  Concise  Introduction  to  the  Theory  (New  York: 
McGraw-Hill  Book  Co.,  1969),  pp.  266-273. 
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Substituting  for  pj(a-q)  in  the  momentum  equation  from  the  continuity  equation 
yields 


Pf  -  Pi 


Ap  ■  p {a2  -  p  £a(a-q  ) 
4p  “Pi*  (a  ••  (a-q)j 
4p  -  P£aq  , 


which  is  the  Rankine-Hugoniot  relation.  If  p  is  taken  to  be  the  gauge  shock 
pressure,  and  since  pi  *  atmospheric,  then  Ap*p.  Dropping  the  subscript  on 
the  initial  or  undisturbed  mass  density,  the  Rankine-Hugoniot  relation  becomes 
simply 


P  "  P«q  , 


(13) 


the  impact  shock  pressure. 


Determination  of  Shock  Wave  Speed  in  Water 

Shock  wave  speed  in  water,  a,  increases  with  the  associated  shock  pressure, 

with 


a+c  as  p  •*  0. 


Values  of  shock  wave  speed  as  a  function  of  shock  pressure  have  been  determined  by 
measurement  on  shocks  in  water.  The  paper  by  Rice  and  Walsh^  tabulates  the 
results  of  this  experimental  Hugoniot  curve  data  for  water  initially  at  20°C  and 
1  atmosphere  to  shock  pressures  of  250  kilobars.  A  plot  of  this  data  to  20 
kilobars  is  given  by  Figure  7. 

Computational  use  of  this  data  table  is  facilitated  by  use  of  a  four  point 
Lagrangian  interpolating  polynomial,  a  description  and  formulation  of  which  is 
given  in  C.F.  Gerald's  text.  If  p  is  a  shock  pressure  which  lies  within  the 
table  as:  P1<P2<P<P3<P4 »  w^ere  Pl»  P2>  P3>  and  P4  Hre  successive  shock  pressure 
table  entries  with  a^,  &2>  ®3i  and  a4  as  their  corresponding  shock  wave  speeds, 
then  the  interpolant  is  given  by  the  cubic 


2Rice,  M.H.  and  Walsh,  J.M.,  "Equation  of  State  of  Water  to  250  Kilobars,"  The 
Journal  of  Chemical  Physics,  Vol.  26,  No.  4,  1957,  p.  824. 

^Gerald,  C.F.,  Applied  Numerical  Analysis  (Reading:  Addison-Wesley  Publishing 
Co. ,  1978) ,  pp.  174-176. 


21 


NSWC  TR  82-438 


IBBBVNM 

■■iuan 


SSSSSSB! 

iSSSSSS 


Sanaa 

■nun 


SSBSSSSSBSSi 

SES5SSSS8SS! 

■■■■■■■■■■■I 


l8BBBBRRBSB8«|iRBBBSHSEBBBBBSBHE»BSSH8BSS&BS8S 

BnBBflBBBBBBBBBifllBBBBiBiflBBBBBEBflBaBBflfliBBBBi5iiBBBaBBiiBni 
■BBBBBBlBBBBaBBBBBgBBEBBBagBaflBBBBBBBBBBBBBBBaBfiBBEBBBBBBI 


gBBBBBBlBBBBBBaBaBBBBgaBaaBBBBBBBaBBBBBBBBBBBBSa3BBB3aaaBaB| 

aaaaBainaBaaBaaBaaiaBaaaaaiaaBaaBaBBaBgawaaaBaaaaaaaaaaBaal 

^^^^■BBflgBaflBBflaBBaBBaflaaBaBflaBBBBBflO|BS0!i«BBaaaMii 
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBCMBBBBiBBBBBWwBai 


bbSbSSSiSSBbSSI 

BBBBBBr  I  I  I  *  ] 

BBBBBBBigiBBBWiwal 

laSSSSaSEBSkSSSaS!  Si  iSBBBS 

IBBBBBBBBBBB IBBBBBBBBBBBBBBBflBBBBBBBBBBBB 

alBBBBBBBBBUBBNBI  II  It  IBBBflBBBBBBK^MBiglBB 
■BNBBNfliNBBrBBNBNBBI-  BBBBB BB  II IBBI  IBBBBBBB 
iBBBBBBBBBBBBiBBflBBBBBBB  I  IBBBB  -  I  ■ 

I  B  •  IB  1BBBBBBB  ■  BflBBI  IBBBB 

liiBBBBBBBBBailBBBBBBBBBBBBBflBBBBBBflflBBBBBl  I 

aUBflBBBBBflflBIBBBBBBBBBBBBBBBflBBBBagBBBBBEXa 
iBBBBBaflBBBBBlBBBBBBBBBBaBBBBBBBBBiBBBBBBBa^^^^^^^^^^^^H 

IbbbBSbbBbS  BBBiBBbBBBbBbBBBbbBbbBBbBbBBbBBBBbbbbbbBbbbbbbbbbbI 

iBaBBBBBBBBBBBBIBaBBBBBBaBBBBBBBBBaBBBBBBBBaBBfiSBaBBBBBiBBBBBl 
iBBifiBBflflBBBBflflllflBflflBiBBBBBBBBBBflBBBBBBBBBBBiiBBBBBBflBBflBiBBBl 
IflflBBflBBBBBBBBBI'BBflBBBBBBBBBflBBBBBBaBflBBiaBBBBBBBBflBBflfliBBiaBl 
lBBBBBBBBBflBBBBB«BBBBflBBiBBBBBflBBBBBBBBBBBflBBBfliBBBBS9lBBBBBE| 
|BBBBBBBBflBBBBBBIilBBBBflBBBBBBBBBBBBBBflBBBBfliBflBBiBBflBBBflBBBBBB| 
IbBBBBBBBBBBbBbBI'BBBBBBBBBB ■■■■■■BBlBaBBgBBBBBBBBBBBBBBBBBBBBl 
|BBBBflBBBflBflBBBBBlBBBBBBBBBBBBBBflBBBBBBaBBBBflBBBBBBBBBflBBBBBBl 


ib&bbbbbbS| 

iBflflBBBBBBBBBBBflBBBBBBBBBBBBBBBBiBBBBBBBBBBBBBBBaBBBBBBBBBfiaBj 

Ibbbbbbbbbbbbbbbbbbbb.ibbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbm 

liBaBBBBBBBBBBBBBBBBBIBBBBBaBBBBBBBBBBBBBaflBBBBiBBBaBBBBBBBBl 

SIwbbbbbbbbbbbbbbbbbbbbibbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbBbBbbbi 
Bbbbbbbbbbbbbbbb^bbb^bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbi 
Hbbbbbbbbbbbbbbbbbbbbbbbwbbbbbwbbbbbbbbbbbbbbbbbbbwbwwbbwbbi 

llBBBBBBBBS  niBBBflflBI  IfllifllBt  iflflifliflBBBBBBBBIBBRBi 


FIGURE  7.  SHOCK  WAVE  SPEED  IN  WATER  AS  A  FUNCTION  OF  SHOCK  PRESSURE 
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(p-p?)  (p-pJ  (p-pJ  (p-p,)  Cp-p« )  (p-p,) 

a(p)  “ - - - J  H  ,  ai  +  1  -  J _ z _ a2 

(P1-P2HP1“P3KP1-P4)  <p2“pl)(p2"p3)‘p2_p4) 


(p-p  )  (p-p.)  (p-p,)  (p-p,)  (P-P,)  <P~P,) 

+  - L - £ - ft —  a3  +  - 1 - 2 - J —  a4. 

(P3-P1)(P3-P2^P3-P4^  (P4-P1^P4"P2)(P4_P3) 

Solution  of  the  Rankine-Hugor.L'.'t  Snock  Relation 

Since  the  shock  wave  sp»ed  is  an  empirical,  function  of  the  shock  pressure, 
and  further,  the  water  particle  speed  is  a  function  of  missile  face  deforntation 
rate,  and  hence  the  shock  pressure,  the  Rankine-Hugoniot  relation, 


P  -  p  a(p)q(p)  , 


(R13) 


■is  an  implicit  function  of  shock  pressure. 

Consider  the  following  function  of  some  trial  pressure  defined  by 

G(pi).  ■  Px  -  pa(p1)q(p1). 


This  is  called  the  shock  pressure  error  function  Since  a  solution  of  Equation 
(R13)  for  the  shock  pressure  is  a  root  of  the  error  function,  or 


G(p)  ■  p  -  paq  -  0. 


(15) 


A  root  is  found  by  an  iteration  of  trial  pressures  until  Equation  (15)  is 
satisfied  to  within  some  specified  tolerance.  One  iterative  technique  used  for 
converging  on  a  root  of  an  implicit  function  is  the  Secant  Method,  which  is 
described  in  a  text  by  Atkinson,^  with  its  first  steps  procedure  illustrated  by 
Figure  8. 

The  initial  trial  value  of  shock  pressure  used  by  the  method  is  the  weak  or 
sonic  shock  pressure  produced  by  a  rigid  body  entering  the  water  with  speed  V'. 
For  this  case  the  Rankine-Hugoniot  relation  reduces  to  the  explicit  form 


Pi  “  pcV' . 


(16) 


^Atkinson,  K.E.,  An  Introduction  to  Numerical  Analysis  (New  York:  John  Wiley  & 
Sons,  1978),  pp.  48-52. 
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Since  c<_a  and  q<V'  this  trial  pressure  is  a  lower  bound  on  the  convergent 
shock  pressure.  A  second  trial  value  is  needed  to  start  the  iteration  process  and 
this  is  supplied  simply  by  an  arbitrary  102  increase  in  the  first  trial  value, 
i.e. , 


p2  -  Pi  +  0.1  Pi. 


(17) 


Between  these  two  values  on  the  G(p)"0  curve  a  secant  is  extended  to  intersect  the 
curve  at  a  point,  determining  the  next  trial  pressure  P3.  This  process  is  then 
repeated  using  the  latter  two  trial  values  to  determine  successive  values 
according  to  the  recursion  formula 


Pi+1 


Pi-1 


0<Pi-i>' 


G(pi)  >  1  ■  3,  4,  5,  n 


(18) 


until  convergence  to  within  a  specified  accuracy  on  the  ratio  change  in  two 
successive  values  occurs,  i.e.  until 


Pn+1  I  Pn 
pn 


<  error  specified 


(19) 


EXPANSION  WAVE  SPEED  ACROSS  DEFORMING  MISSILE  FACE 

During  water  entry  the  shock  pressure  behind  the  compression  shock  wave  front 
is  then  later  relieved  back  to  the  initial  pressure  with  the  passage  of  a 
rarefaction  or  expansion  wave  generated  by  the  elastic  recovery  of  the  water. 

This  expansion  wave  however  travels  only  at  the  sonic  or  weak  disturbance  speed  c' 
since  water  cannot  sustain  a  tensile  stress.  Moving  outward  from  the  initial 
contact  point  this  expansion  wave  moves  across  the  deforming  face  of  the  missile 
with  a  component  speed  b<c ' ,  since  this  face  is  not  parallel  to  the  local  ray  of 
the  expansion  wave.  This  component  of  expansion  wave  speed  is  determined  by  the 
trigonometry  of  Figure  9  as  the  face  side  bt  opposite  the  expansion  wave  ray  c't 
and  the  extended  deforming  entry  side  V't-c2emt. 

The  angle  6  between  the  extended  deforming  face  bt  and  the  entry  side 
V't-tjja^t  is  the  supplement  of  the  complement  of  the  deformation  angle  u, 
as  can  be  seen  from  the  triangles  in  Figure  9,  or 


6  *  it  -  (JL  -  hi  )  , 

2 

6  *  I.  +  w  •  (20) 

2 
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The  angle  41  between  the  extended  deforming  face  bt  and  the  expansion  wave  ray 
c't  is  found  by  use  of  the  Sine  Law  as 


Sin-1 


(V'  -  e  a  )  tin  6 
z  m 

- - 


(21) 


This  leaves  the  angle  f  between  the  expansion  wave  ray  and  entry  side  as  the 
triangle's  supplement,  or 

Y  ■  it  -(«+♦).  (22) 

Finally  employing  the  Sine  Law  again  yields  the  deforming  face  component  of 
expansion  wave  speed  as 

b  -  Si  8in  ♦  .  . 
sin  6 


The  projection  upon  the  rigid  missile  face  of  this  speed  is 


c  M  b  cos  ui. 


Thus  the  component  of  expansion  wave  speed  moving  across  the  undeformed 
circumferential  circle  of  the  missile's  rigid  face  is 


c  ■  c  1  s  in  <p  cos  u 
sin  £ 


(23) 


This  is  the  speed  required  for  determining  impact  Bhock  duration  time,  expansion 
wave  profiles,  and  shock  pressure  areas  relative  to  the  missile's  initial 
dimensions  and  geometric  configuration,  yielding  results  which  will  be  projected 
axial  components. 

Note  that  for  the  special  case  of  normal  entry  impact  (a“0  implies 
was  previously  shown)  Equations  (20),  (21),  (22),  and  (23)  reduce  to: 
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•*/  /v'  -  c  a  \2 

c  -  o'  ■yi  -  - ,  z  m) 

c  -  ^'2  -  (V1  -  c8am)2  . 

which  is  again  consistent  with  the  results  of  D.  Laumbach's  paper  for  normal  entry 
impact . 


TIME  DURATION  OF  IMPACT  SHOCK 

When  the  expansion  wave  has  engulfed  the  shock  pressure  region  over  the 
entire  face  of  the  missile  the  shock  phenomena  has  ceased  and  only  hydrodynamic 
forces  are  left.  With  the  missile  entering  the  water  obliquely  the  wetting  line 
moves  across  the  rigid  face  at  a  finite  speed  V,  e.g.  see  Figure,  1,  where 


V 


V1 

tan  a 


(24) 


and  with  V  »  as  a  +  0. 
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At  this  speed  the  wetting  line  will  have  traveled  the  missile's  diameter  in  the 
t  ime 


Te 


2r  „  2r  tan  a 
V  V’ 


(25) 


the  water  entry  or  full  wetting  time.  At  this  point  the  initial  expansion  wave 
has  traveled  the  distance  c'fe  leaving  the  remaining  diametral  length  2r  -  cTe 
under  shock  pressure.  This  remaining  length  is  then  engulfed  by  both  the  initial 
and  now  the  final  expansion  waves  coming  in  from  opposite  diametral  edges  in  the 
plane  of  symmetry  as  shown  by  Figure  10.  This  remaining  length  will  be 
transversed  by  the  two  waves  in  the  time  (2r  -  cTe)/2c.  Thus  the  impact  shock 
duration  time  will  be  the  sum  of  this  plus  the  water  entry  time)  or 


Te 


2r  -  cT 
_ e 

2c 


Substituting  for  Te  from  Equation  (25)  and  simplifying  yields 


a  r  tan  a  +  r 
1  V'  ’  c 


(26) 


where  T£  (a"0)  ■  r/c  for  normal  impact,  in  which  case  expansion  waves  are 
generated  on  the  full  circumference  simultaneously  at  t"0  and  there  is  no  wetting 
line  motion. 


SHOCK  AREA  ON  MISSILE  FACE 

In  order  to  determine,  the  impact  shock  force  and  impulse  acting  on  the 
missile  it  is  now  necessary  to  determine  the  shock  pressure  area  over  the  face  of 
the  missile.  This  area  will  be  encompassed  by  the  expansion  wave  profile  on  the 
face  and  the  wetting  line  where  the  shock  region  hsd  its  origin,  when  this 
exists.  This  shock  area  profile  is  now  determined  for  oblique  and  normal  water 
entry.  Note  that  these  profiles  will  be  constructed  on  the  rigid  face  as  a 
projection  of  the  deformed  face  with  the  resulting  shock  areas,  and  hence  forces, 
being  axial  components  of  the  deformed  face  normal. 


Shock  Area  Profile  for  Oblique  Entry 

For  oblique  entry  a  wetting  line  is  established  which  moves  across  the 
undeformed  face  of  the  missile  at  finite  speed  V  given  by  Equation  (24),  Chasing 
this  wetting  line  is  the  expansion  wave  profile  moving  at  speed  c,  given  by 
Equation  (23),  with  expansion  wave  rays  being  continuously  generated  on  the 
circumference  at  the  cusps  of  the  wetting  line.  When  the  wetting  line  has 
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1 


traveled  2r  to  the  end  of  the  face,  the  face  has  become  fully  wetted  and  full 
water  entry  has  occurred  in  the  time  Te  given  by  Equation  (25).  After  this  time 
no  additional  expansion  wave  rays  are  generated  and  the  remaining  shock  area  is 
engulfed  by  the  then  existing  full  circumferential  expansion  wave  by  impact  shock 
duration  time  T^,  after  initial  water  contact,  given  by  Equation  (26). 

During  Wetting  Entry.  For  an  instantaneous  time  since  initial  water  contact 
of  t<Te  Che  shock  area  envelope  is  bounded  by  a  wetting  line,  as  shown  in 
Figure  11,  as  well  as  the  expansion  wave  profile.  This  profile  ia  an  envelope  of 
the  parametric  family  of  expansion  wave  rays  generated  along  the  circumference 
between  the  initial  contact  point  and  the  current  cusps  of  the  wetting  line.  From 
the  initial  contact  point  the  ray  which  was  generated  at  parametric  time  t“0, 
now  has  the  instantaneous  radius  p(t*0)  *  ct.  At  the  wetting  line  circumferential 
cusp  an  expansion  wave  ray  ia  just  being  generated,  thus  ita  radius  is  p(T»t)  ■  0. 
So  that  at  an  intermediate  parametric  time  t,  CKr^t,  the  wetting  line  occupied 
the  position  (£,n)  on  the  circumference,  having  traveled  the  distance  Vt , 
and  generated  the  expansion  wave  ray  whose  radius  at  current  time  t  ia  now 


p  (t  )  ■  c(t-t )  (27) 

I 

with  a  missile  circumferential  origin  in  the  (x,y)  face  coordinates  of 

a 

5  ■  r  -  Vt  and  (28) 

n  M  Vr^  -  £  2  .  (29 ) 


These  expansion  wave  rays  generate  the  parametric  family  of  circles! 
(x  -  ?  )2  +  (y  -  n)2  -  P2  ■  0  , 


The  expansion  wave  profile  on  the  ri{jid  missile  face  will  then  be  the  curve 
tangent  to  this  family  of  expanding  circles.  The  necessary  conditions  which  such 
an  envelope  E,  cf  a  family  of  curves  C:  f(x,y,T)  ■  0,  must  satisfy  are 
developed  in  R.  Courant's  text.  0  The  envelope  must  satisfy  the  system  E: 

1)  f (x,y ,t  )  *  0  and 

2)  fT  (x,y  ,t  )  *  0  . 


l^Courant,  R. ,  Differential,  and  Integral  Calculus,  Vol.  II  (New  York: 
Interscience  Publishers,  1936),  pp.  171—179. 
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A  V 


FIGURE  117  SHOCft  AREA  ENVELOPE  DURING  WETTING  ENTRY 
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with 

f(x,y,i)  ■  (x  -  5)2  +  (y  -  n)2  -  p2  *  0 

then 

il  ■  2(x  -  t)  *£)  +  2 (y  -  n)  (-  51)  -  2p  ®£  -  0, 

31  3i  3t  3t 


where 


2£ 

3t 


-  -V 


» 


_  J, 

iH  ■  I  (r2  -  £2)  2  .  (-2)C  ii  or 

3t  2  3t 


®1  -  I  V ,  and 
3t  n 


3  P 
3l 


-c . 


Thus 

ii  ■  (x  -  O  +  (y  -  n)  (-  £•)  -  -  o 

3i  n  v 


or,  with  4 


<30; 


the  necessary  conditions  for  the  envelope  E  become: 

1)  (x-£  )2  +  (y-n)2  -  p2  ■  0  and 

2)  x-JLy+4  "0  . 

n 
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Solving  condition  (2)  for  y  gives 


y  (t  ) 


H  (x  +  5 )  . 
? 


(31 


Then  substituting  this  expresion  for  y  into  condition  (l),  collecting  like  terms, 
and  making  use  of  the  relation  52  +  n2  “  r2  yields  the  quadratic  in  x: 

r2x2  -  2(r25  -  n2«  )x  +  [r2r, 2  -  n25  (25  -  S  )  -  52«2]  "0. 


With 


A  «  r2, 

B  ■  -2(r25  -  n2i ) ,  and 
C  -  r25  2  -  n2«  (25  -  6)  -  5  2P2 
the  roots  will  be 

"•2 '  -  {k) -J{k)2 '  i  ■ 


where 


2  2 

£-iu2UL-i.  and 
r2 


(L.)2  -  £  -  *2<*V  ~  nV)  >  o. 

\ 2A/  A  r4 


This  discriminate  will  be  positive,  and  hence  x(i  )  real,  if  p>6  (where  r>n  by 
definition)  or,  from  Equations  (27)  and  (30),  if  V>c.  This  is  the  wetting  line 
speed  condition  required  for  a  nonzero  shock  pressure  area  to  be  established,  cf. 
Equation  (9). 
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Of  the  two  real  roots,  x^  and  x^ ,  the  smaller  is  the  one  of  interest 
since  it  lies  within  the  missile's  circumference.  Thus 


x(x)  -  C  "  3_*_ 


*  £  Vr 


n — rr 

o  -  n  £ 


(32) 


Equations  (32)  and  (31)  are  the  parametric  coordinates  of  the  expansion  wave 
profile  projected  over  the  missile's  rigid  face,  [x(t),  y  Ct  ) }  for  0<T<t,  at 
a  current  instantaneous  impact  time  t  since  initial  water  contact. 

After  Face  is  Fully  Wetted.  For  t>Te  the  missile's  face  is  fully  wetted, 
the  wetting  line  no  longer  exists,  and  established  expansion  waves  emanate  from 
the  full  cii*cumference  of  the  face  as  shown  by  Figure  12.  Since  no  n* v  expansion 
wave  rays  are  generated  the  established  ones  will  have  their  circumferential 
origins  as  previously  described  by  Equations  (28)  and  (29)  when  t**Te,  or 


U(t),  n(t)]  ,  0  <r<Tfi  for  t>Te  . 


With  the  waves  continuing  to  expand  from  their  established  positions  with  speed  c, 
Equation  (27)  for  their  parametric  radius  remains  unchanged  as 


p(x>_c(t-x)  ,  0<x£Te  for  t>Tfi . 

Hence  the  expansion  wave  profile  [x(T),y(x)]  is  described  by  the  same  set  of 
equations,  Equations  (32)  and  (31),  as  for  the  wetting  entry  but  with  the  time 
parameter  limited  accordingly  as 


for  t<Te 
for  t>Te 


(33) 


Integration  for  the  Shock  Pressure  Area 

With  the  envelope  of  the  shock  pressure  area  on  the  missile's  rigid  face — 
the  expansion  wave  profile  described,  this  region  may  now  be  integrated  to 
determine  the  projected  shock  pressure  area  on  the  rigid  face.  This  integration 
is  performed  numerically  with  equal  time  steps  in  the  time  parameter  x  describing 
the  profile.  These  time  steps  are  a  sufficiently  large  division  of  the  current 
impact  time  t,  or  the  water  entry  time  Te  if  t>Te,  so  as  to  insure 
sufficient  numerical  accuracy.  Since  equal  steps  in  the  parameter  of  a  parametric 
curve  do  not,  in  general,  result  in  equal  divisions  of  the  x-ordinate,  the 
standard  numerical  integration  schemes  based  upon  equal  interval  divisions  of  the 
ordinate  cannot  be  used.  Integration  is  instead  based  on  a  generalization  of 
Simpson's  Rule  for  parametric  curves,  illustrated  by  Figure  13,  in  which  a 
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Lagrangian  second  order  interpolating  polynomial  P2(x)  is  passed  through  three 
successive  points  on  the  parametric  curve,  approximating  that  curve  and 
corresponding  to  a  single  integration  panel. 

If  (x^.y^),  (x2,y2) ,  and  (xj^)  are  three  successive  points  on  the  expansion 
wave  profile  curve  corresponding  to  the  successive  time  parameters  tj,  t2  “  t  j.  + 

A-r ,  and  13  ■  t  ]_  +  2At  respectively,  then  the  Lagrangian  interpolating  polynomial 
for  these  three  points  is 


P2(x) 


(x-x„)  (x-x, )  (x-x,  )  (x-x,)  (x-x,  )  (x-x, ) 

2  J  y  ■  I  4.  1  ^  y «  H*  ^  *•  y  ^ 

Txr-x^TTxp-xjy  (x2--x^)  (x2-X3)  (X3-X3)  CX3-X2) 


Integrating  this  for  the  upper  half  profile  curve,  with  x(t)  and  yOt)  provided 
by  Equations  (32)  and  (31)  reopectively ,  for  Ckr<t,  or  Te  if  t>Te,  yields 
and  approximation  of  half  the  shock  pressure  area,  because  of  symmetry  about  the 
x-axis,  for  a  single  panel,  or 


P2(x)dx. 


which,  after  carrying  out  the  integration,  results  in  the  area  approximating 
formula 


fl  «  3  ^x3  "  xl^  “  1  ^x2  f  x3^x3  ~  xl^  *  X2X3^X3  "  xl^ 

2~  (hi  -  xjKxI  -  X3) 


yi 


•J  (xlj  -  xj[)  -  (x^  +  x3)(x|  -  x2)  +  x^XjCx-j  - 

(x2  X1)(x2  -  X3I 


?2 


<x3  - 


x^ )  -  (x^  + 


x2)(x2 


"(X3  -  x^)(x3  -  x 2 ) 


x2) 


xlx2(x3  ~ 


Xl) 


y3 


(34) 


Since  two  increments  of  the  time  parameter  are  required  per  panel  the  number  of 
division  steps  must  be  even.  Thus  for  a  given  instantaneous  impact  time  t 
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S  . 

2 


n 


AS. 

~  »  0<«ti  . 


i"l»3 ,5, . . . 

where  i  ■  1  corresponds  to  t  *  0 

.  i 

i  *  n  corresponds  to  i  ' 

with  At  m  t/n  or  Te/n  and  n  is  odd 
intervals. 


and 

t  for  t<Te 
Te  for  t>Te 

i.e.  there  are  n-L  even  time  parameter 


The  Special  Case  of  Normal  Entry 

For  the  case  of  normal  water  entry  impact,  a"0,  expansion  waves  are 
generated  on  the  full  circumference  of  the  missile's  face  instantly,  at  t*0,  and 
simultaneously.  This  results  in  a  concentric  circular  expansion  wave  profile 
shrinking  the  shock  pressure  area  with  speed  c.  Thus  the  shock  pressure  area  for 
this  case  is  simply  the  area  of  a  circle  with  a  radius  diminishing  at  the  rate  c, 
or 


S(a-O)  “it  (r-ct)2  ,  0<t<Ti  (35) 

IMPACT  SHOCK.  FORCE  AND  IMPULSE 

Finally  with  the  impact  shock  pressure  and  instantaneous  projected  shock 
pressure  area  determined,  the  axial  component  of  the  impact  shock  force  acting  on 
the  missile  is 

F(t)  *  pS(t) .  (36) 

With  this  result  the  total  axial  shock  impulse  may  also  be  found  as 


I 


/  F(t/dt. 
o 


This  integral  is  evaluated  numerically  on  equal  intervals  of  t  by  use  of  Simpson's 
1/3  Rule: 
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n 

1  -  *-  ^7  CFi  ■+  4Fi+l  +  Fi+2)  ,  (37) 

r*l,3,5,... 


where  i  K  1  corresponds  to  t  “0, 
i  -  n  corresponds  to  t  *  T^, 
and  At  “  Tj/n  tor  n  odd. 

As  before  n-1  even  time  intervals  are  required  for  the  approximate  integration, 
with  a  larger  value  of  n  improving  the  accuracy  on  the  total  impulse. 

IMPACT  SHOCK  DYNAMICS  COMPUTATION  PROGRAM  ’IMPACT' 


Computation  of  the  shock  dynamics  for  a  flat  faced  cylindrical  missile 
impacting  the  water's  aurfsce  is  carried  out  using  the  numbered  equations 
developed  in  the  previous  section  by  the  program  IMPACT.  A  listing  of  this 
program  is  given  in  Appendix  A.  The  programming  language  used  is  CDC's  version 
of  FORTRAN  IV.  The  program  is  then  applied  to  a  hypothetical,  but  hopefully 
typical,  example  problem  of  a  missile  entering  the  water  at  various  angles. 
Output  for  this  example  is  given  by  Appendix  B.  The  results  of  this  example  are 
then  discussed  with  regard  to  water  impact  physics  and  the  influence  of  entry 
angle.  , 


DESCRIPTION  OF  PROGRAM  STRUCTURE 

The  program  IMPACT  has  been  written  in  a  simplistic  or  "engineering"  format 
with  emphasis  upon  clarity  and  ease  of  modification  rather  than  efficiency.  The 
program  structure  is  linear,  following  the  same  sequence  of  steps  used  in  the 
derivation  of  the  basic  equations  for  impact  shock  dynamics*  An  outline  of  the 
program  code  is  given  below,  with  the  program  symbolic  variable,  its  correspond¬ 
ing  equation  symbol,  and  its  equation  number  in  the  previous  sections  included. 

I.  CONSTANT  PHYSICAL  DATA  ENTRY  (Automatic) 

A.  Rice  &  Walsh  Shock  Wave  Speed  Data 

1.  Shock  pressure,  D(1,I)  *  p,  psi 

2.  Shock  wave  speed,  D(2,l)  ■  a,  in/sec. 

B.  Undisturbed  Properties  of  Water 

1.  Geometric  pi,  PI  ■  it 

2.  Density,  RHO  ■  p,  lb-secI. 2 3/in^ 

3.  Speed  of  sound,  Cl  ■  c ' ,  in/sec. 
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II.  ENTRY  OF  INDEPENDENT  VARIABLES  (Prompted) 

A.  Missile  Radius,  R  ■  r,  in 

B.  Impact  Velocity,  V2  ■  V1,  fps 

C.  Entry  Angle,  ALPHA1  *  a,  deg. 

D.  Missile  Weight  Density,  RH0M1  ■  Pm/g»  lb/in"* 

E.  Young's  Elastic  Modulus,  E  ■  E,  ps i 

F.  Poission's  Ratio,  NU  ■  v 

G.  Number  of  Time  Intervals  for  Printout,  M  *■  n/2 

III.  PRINTING  OF  TABLE  HEADER 

A.  Title 

B.  Input  Data  Review 

IV.  COMPUTATION  OF  CONSTANT  SHOCK  DYNAMICS  RESULTS 

A.  Missile  Elastic  Wave  Speed,  AM  “  ara  (3) 

B.  Shock  Pressure  Computation 

1.  Initial  Trial  Shock  Pressure,  Pi  -  pi  (16) 

2.  Subroutine  PRESS: 

a.  Axial  missile  strain,  EZ  **  t  z  (4) 

b.  Missile  face  deformation  Angle  OMEGA  “id  (5) 

c.  Resultant  water  particle  speed,  U  ■  u  (6) 

d.  Subroutine  ITABLE:  Interpolation  for  water  shock  wave 
speed  data  from  trial  shock  pressure  (14) 

e.  Check  on  surface  disturbance  speed  for  shock  (8) 

f.  Mach  angle,  BETA  ■  13  (7) 

g.  Normal  water  particle  speed,  Q  ■  q  (12) 

h.  Rankine-Hugoniot  shock  pressure,  PS  “  p  (13) 

i.  Shock  pressure  error  function,  G  ■  G(p)  (15) 

3.  Second  Trial  Shock  Pressure,  P2  ■  p2  (17) 

C.  Subroutine  SECANT:  New  iterate  value  for  shock  pressure  (18) 

D.  Expansion  Wave  Speed  Across  Rigid  Face,  c 

1.  DELTA  -  6  (20) 

2.  PHI  -  *  (21) 

3.  PSI  -  «>  (22) 

4.  C  -  c  (23) 

E.  Time  at  Full  Wetting  Entry,  TE  “  Te  (25) 

F.  Time  Duration  of  Impact  Shock,  TI  ■  Tj  (26) 

V.  PRINTOUT  OF  CONSTANT  SHOCK  DYNAMICS  RESULTS 

A.  Shock  Wave  Speed  in  Missile,  AM  *  an,  in/sec. 

B.  Elastic  Strain  on  Missile  Face,  EZ  "  t z,  in/in 

C.  Deformation  Angle  of  Missile  Face,  OMEGA!  *  w,  deg. 

D.  Water  Particle  Velocity  in  Shock,  U  ■  u,  in/aec. 

E.  Shock  Wave  Speed  in  Water,  A  ■  a,  in/sec. 

F.  Mach  Angle  of  Shock  Wave  in  Water,  BETA1  ■  13,  deg. 

G.  Water  Particle  Speed  Normal  to  Shock  Wave,  Q  *  q,  in/Bec. 

H.  Rankine-Hugoniot  Shock  Presoure,  P  ■  p,  psi 

I.  Expansion  Wave  Speed  Across  Missile  Face,  C  ■  c,  in/sec . 

J.  Time  at  Full  Wetting  Entry,  TE  ■  te,  sec. 

K.  Time  Duration  of  Impact  Shock,  TI  ■  T^,  sec. 
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VI.  COMPUTATION  OF  TIME  DEPENDENT  RESULTS 

A.  Time  Step  Interval,  H  ■  At 

B.  Initialize  Impulse,  FT  ■  I  ■  0. 

C.  Shock  Pressure  Area  Determination,  S(t) 

1.  For  Normal  Impact,  S  *•  S  (35) 

2.  For  Oblique  Impact: 

a.  Wetting  line  speed,  V  ■  V  (24) 

b.  Check  on  wetting  line  speed  for  shock  (9) 

c.  Subroutine  AREA:  Half  area  between  wetting  line  and 
expansion  wave  profile 

i.  Initialize  area,  A  ■  S/2  ■  0. 

ii.  Determine  time  parameter  limit,  T1  (33) 

iii.  Time  parameter  for  profile  curve,  TAU  ■  t 

iv.  XI  -  5  <28) 

v.  ETA  *  n  (29) 

vi.  RHO  -  p  (27) 

vii.  DELTA  -  4  (30) 

viii.  Coordinates  of  expansion  wave  profile  curve: 

X  «  x(t )  (32) 

Y  -  y(t)  (31) 

ix.  Simpson's  Rule  integration  for  unequal  ordinate 

intervals,  A  ■  S/2  (34) 

D.  Impact  Shock  Force,  F  “  F(t)  (36) 

E.  Subroutine  SIMP:  Integration  of  impact  shock  force  over  time  for 
impulse,  FT  ■  I(t)  (37) 

VII.  PRINTOUT  OF  SHOCK  DYNAMICS  RESULTS 

A.  Water  Entry  Time,  T  **  t,  sec. 

B.  Shock  Pressure  Area,  S  “  S,  in2 

C.  Impact  Shock  Force,  F  ■  F,  lb 

D.  Impact  Shock  Impulse,  FT  ■  I,  lb-sec. 

The  error  tolerance  on  convergence  to  a  root,  used  with  the  subroutine 
SECANT,  as  a  ratio  change  in  its  previous  value  is  0.0001.  For  the  numerical 
integration  of  shock  pressure  area,  carried  out  in  the  subroutine  AREA,  200 
subdivisions  of  the  parametric  time  interval  is  used.  This  has  been  found  by 
trial  to  result  in  a  lower  bound  on  area  with  an  error  of  about  5  parts  in 
100,000.  Finally,  the  number  of  steps  used  for  integration  of  the  total  shock 
impulse  is  twice  the  number  of  time  intervals  of  shock  duration  time  requested 
of  the  user  for  the  printing  of  the  time  dependent  results.  With  this  integra¬ 
tion  carried  out  by  Simpson's  Rule  in  the  subroutine  SIMP,  the  error  will  be, 
from  Gerald,  C.F.,  "Applied  Numerical  Analysis",  p.  215, 


t.  /t.v*  . 

Error  -  -  Fiv(t)  ,  0<t<Ti 

T50  \ 2M/  —  1 
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where:  is  the  impact  duration  time, 

M  the  number  of  time  intervals  for  printout  requested  by  the  user,  and 

F1VU:)  is  the  fourth  time  derivative  of  the  impact  shock  force  (which 
unfortunately  is  not  known  prior  to  computation)* 

The  program  makes  checks  on  the  conditions  necessary  for  the  impact  to  be 
shock  producing  in  accordance  with  Equations  (8)  and  (9).  If  either  of  these 
conditions  is  not  satisfied  then  the  appropriate  message  is  printed  and  execution 
halts*  Thus  the  progam  IMPACT  will  give  shock  dynamics  results  only  for  those 
entry  angles  which  are  less  than  the  critical  entry  angle. 


USE  OP  THE  PROGRAM 

The  program  has  been  written  for  interactive  use  on  a  time  sharing  terminal 
with  the  program  normally  kept  in  the  users  library  as  an  active  file,  or  loaded 
into  it  by  a  card  deck.  All  the  necessary  independent  variables  needed  to 
describe  a  particular  water  entry  impact  problem  are  prompted  for  by  the  running 
program  in  typical  engineering  units,  specified  in  the  prompt  as: 

o  ENTER  RADIUS  OF  MISSILE,  inches 

o  IMPACT  VELOCITY  OF  MISSILE,  feet  per  second 

o  ANGLE  OF  ENTRY  FROM  VERTICAL,  degrees 

o  ENTER  WEIGHT  DENSITY  OF  MISSILE,  pounds  per  cubic  inch 

o  YOUNG'S  ELASTIC  MODULUS ,  pounds  per  square  inch 

o  POISSION'S  RATIO  FOR  MISSILE,  (dimensionless) 

In  addition  to  the  independent  variables  the  number  of  time  intervals  in 
the  impact  shock  duration  time  for  print  out  is  requested  of  the  user  as: 

o  ENTER  THE  NUMBER  OF  TIME  INTERVALS  DESIRED  - 

Since  the  integration  of  shock  force  over  time  is  accomplished  by  the  numerical 
method  of  Simpson's  Rule,  a  larger  number  of  intervals  entered  here  yields  a 
better  computational  accuracy  on  the  total  shock  impulse.  With  the  integration 
steps  used  internally  being  twice  the  number  of  printing  steps  prompted,  an 
entry  of  about  SO  time  intervals  is  generally  sufficient  for  engineering 
accuracy.  An  increase  to  100  resulted  in  only  a  0.01Z  change  in  total  impulse 
for  the  example  with  a  ■  3°. 

With  the  number  of  time  intervals  for  printout  entered  the  program  proceeds 
without  further  action  by  the  user.  The  input  data  is  now  tabulated  and  the 
constant  shock  dynamics  parameters  are  computed,  including  the  Rankine-Hugoniot 
shock  pressure.  During  this  time  a  check  is  made  on  the  surface  disturbance 
speed  as  a  condition  for  an  impact  shock  to  exist.  This  must  be  greater  than 
the  shock  wave  speed  in  accordance  with  Equation  (8).  If  this  is  not  satisfied 
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a  message  is  printed  indicating  so  and  program  execution  halts.  Otherwise  the 
program  proceeds  with  a  printing  of  the  constant  shock  dynamics  results  and 
begins  computation  for  the  time  dependent  results.  Prior  to  determination  of 
the  shock  pressure  area  a  check  is  also  made  on  the  wetting  line  speed.  This 
must  be  greater  than  the  expansion  wave  profile  speed  for  a  shock  pressure  area 
to  exist  as  specified  by  Equation  (9).  If  this  condition  is  not  met  a  message 
is  printed  and  the  program  stops.  Otherwise  execution  proceeds  with  computation 
and  printing  of  the  time  dependent  results  for  each  increment  of  the  impact 
shock  duration  time:  shock  pressure  area,  force,  and  impulse;  with  the  total 
shock  impulse  being  the  last  value  of  impulse  printed. 

A  sample  terminal  session  in  the  use  of  the  program  IMPACT  is  given  by 
Appendix  B  for  the  example  problem  with  a  5°  entry  angle.  Included  are  the 
prompts  for  input,  the  input  data  used,  and  the  resulting  output  tables  of 
constant  and  time  dependent  impact  shock  dynamics. 


EXAMPLE  PROBLEM 

For  the  purposes  of  illustrating  the  use  of  the  program  and  also  for 
verification,  by  the  physical  response,  the  assumptions  used  in  developing  the 
basic  equations,  a  hypothetical  impact  problem  is  solved  with  IMPACT.  This 
example  problem  is  the  water  entry  impact  of  a  flat  faced  cylindrical  missile  of 
aluminum  with  a  6  inch  radius  entering  the  water  at  a  speed  of  1,300  feet  per 
second.  The  program  IMPACT  will  then  be  run  for  various  entry  angles  so  that 
shock  force  and  impulse  as  functions  of  entry  angle  as  well  as  of  time  can  be 
found.  The  necessary  input  data  required  for  this  problem  will  then  be: 

Missile  Radius  «  6.0  in, 

Impact  Velocity  ■  1500  fps , 

Weight  Density  of  Missile  *  0.1  lb/in-*, 

Young's  Modulus  ■  10E6  pai,  and 
Poission's  Ratio  *  0.333; 

with  the  entry  angle  varied  for  each  run  from  0  to  12  degrees. 

RESULTS  OF  THE  EXAMPLE 

The  impact  shock  force  as  a  function  of  water  entry  time  has  been  plotted 
for  entry  angles  of  0,  1,  2,  5,  and  10  degrees  in  Figure  14.  From  this  graph  it 
is  seen  that  maximum  shock  force  occurs  for  normal  entry  (0  deg.  entry  angle)  at 
the  instant  of  water  contact  and  diminishes,  with  an  increasing  time  delay,  with 
increasing  entry  angles.  So  that  as  the  entry  angle  approaches  zero  the  impact 
shock  force-time  profile  approaches  the  normal  entry  case.  The  maximum  impact 
shock  force  is  plotted  for  entry  angle  in  Figure  15.  This  again  illustrates  the 
uniform  decrease  in  shock  force  with  sn  increase  in  entry  angle  from  a  maximum 
for  normal  impact  to,  what  would  be,  zero  for  the  critical  entry  angle,  greater 
than  the  12  degrees  shown.  Figure  16  graphs  the  total  impact  shock  impulse  as  a 
function  of  entry  angle.  This  shows  the  total  impulse  diminishing  with  entry 
angle  as  expected,  though  gradually  at  first  because  of  the  corresponding 
increase  in  shock  duration  time. 


•.  A.’ 'SJiifc 
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ENTRY  ANGLE 


WATER  ENTRY  TIME.  MICRO-SEC 

FIGURE  14.  IMPACT  SHOCK  FORCE  AS  A  FUNCTION  OF  WATER  ENTRY  TIME 


MAXIMUM  IMPACT  SHOCK  FORCE.  10 
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These  results  illustrate  a  physical  response  one  would  expect  of  water 
entry  shock  producing  impacts.  The  results  also  converge  to  those  for  normal 
entry  impact  as  given  by  D.D.  Laumbach's  original  paper.  However,  in  order  to 
determine  whether  the  program  provides  results  which  are  adequately  approximate 
it  is  necessary  to  compare  them  with  experimental  testa  of  oblique  water  impact 
shocks  on  model  missiles. 


POSSIBLE  MODIFICATIONS  AND  IMPROVEMENTS  TO  THE  PRESENT  FORMULATION 


The  most  important  continuation  of  this  study  is  an  experimental 
verification,  for  oblique  entry  shock,  of  the  formulated  model  and  program 
IMPACT.  Beyond  this  the  following  improvements  are  possible  extensions  which 
can  be  carried  out  in  the  simplified  approach  used  thus  far: 

1)  Let  the  entry  velocity  be  variable.  Its  change  is  proportional  to  the 
instantaneous  impulse  being  computed  and  inversely  to  the  missile's  mass. 

2)  Let  the  entry  angle  be  variable  during  impact  shock.  Its  acceleration 
is  proportional  to  the  torque  produced  by  the  Rhock  force.  This  requires 
finding  the  centroid  of  the  shock  pressure  area  and  adding  a  moment  of  inertia 
for  the  missile  to  the  input  data. 

3)  The  assumptions  of  pure  axial  deformation  and  constant  face  deformation 
speed  under  shock  pressure  requires  further  consideration  with  regard  to  the 
theory  of  elasticity  and  wave  mechanics. 

4)  An  analytical  expression  for  the  empirical  shock  wave  speed-shock 
pressure  data  for  water  would  be  a  convenient  substitute  for  the  interpolated 

table . 

5)  A  program  for  determining  the  critical  entry  angle  for  a  given  missile 
entering  the  water  at  a  given  speed  would  be  useful  in  establishing  limits  on 
the  entry  angles  for  shock  producing  impacts. 
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NOMENCLATURE 

a  Shock  wave  speed  in  water  (in/sec.) 

b  Expansion  wave  speed  across  deformed  missile  face  (in/sec.) 

c  Expansion  wave  speed  across  deformed  missile  face  relative,  or 

projected,  to  the  rigid  face  (in/sec,) 

c1  Expansion  wave  speed  in  water,  speed  of  sound  (in/sec.) 

E  Young's  elastic  modulus  of  missile  material  (psi) 

F(t)  Impact  shock  force  on  missile  face  (lb) 

G(p)  Shock  pressure  error  function  (psi) 

I(t)  Impact  shock  impulse  on  missile  (lb-sec.) 

n  Number  of  time  steps  of  impact  shock  duration  time  used  for 

integration  of  impulse 

p  Impact  shock  gauge  pressure  (psi) 

Initial  or  undisturbed  water  pressure  (psi) 

Pi  Final  or  shocked  water  pressure  (psi) 

q  Component  of  shocked  water  particle  velocity  normal  to  water  shock 

wave  front  (in/sec.) 

r  Radius  of  cylindrical  missile  (in) 

S  Shock  pressure  area  on  missile  face  (in^) 

t  Time  since  missile  first  touches  the  water  (sec.) 

Te  Time  at  full  wetting  of  missile's  face  (sec.) 

T£  Time  duration  of  impact  shock  (sec.) 

u  Resultant  shocked  water  particle  velocity  relative  to  undisturbed 

water  (in/sec.) 
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NOMENCLATURE  (Cont.) 

V  Wetting  line  speed  across  missile's  face  in  rigid  face  plane  (in/sec.) 

V'  Missile's  axial  water  entry  velocity  (in/aec.) 

w  Axial  displacement  of  deformed  face  projected  to  initial  entry  side 

(in) 

w  Axial  -speed  of  deforming  face  (in/sec.) 

x(t )  Ordinate  of  expansion  wave  profile  curve  on  rigid  missile  face  (in) 

y(t)  Abscissa  of  expansion  wave  profile  curve  on  rigid  missile  face  (in) 

a  Entry  angle  of  missile's  axis  from  a  normal  to  the  water's  surface 

(rad . ) 

o*  Critical  entry  angle  (rad.) 

13  Mach  angle  of  water  shock  wave  front  (rad.) 

Y  Angle  between  water  shock  wave  normal  and  water  particle  velocity 

vector  (rad.) 

5  Angle  between  deformed  missile  face  and  side  of  missile  in  plane  of 

symmetry  (rad.) 

5(t)  Expansion  wave  profile  curve  translation  parameter  (in) 

dp  Change  in  water  pressure  across  a  shock  wave  front  (lb/in2) 

dp  Change  in  water  density  across  a  shock  wave  front  ( lb-sec^in^) 

e  ,e  Lateral  and  transverse  strains  in  shocked  miscile  (in/ in) 

x  y 

ez  Axial  strain  in  missile  due  to  impact  shock  (in/in) 

n (t )  Ordinate  of  an  expansion  wave  source  on  circumference  of  missile  (in) 

v  Poission's  ratio  elastic  modulus  for  missile 

5 (t )  Abscissa  of  an  expansion  wave  source  on  circumference  of  missile  (in) 

p,Pi  Undisturbed  or  initial  mass  density  of  water  (lb-sec2/in^) 

p (t )  Radius  of  an  expansion  wave  from  an  edge  source  on  missile  face  (in) 

Pf  Final  or  shocked  water  mass  density  ( lb-sec^/in^) 

pm  Undisturbed  mass  density  of  missile  ( lb-sec2/ in^) 

a  ,o  Lateral  and  transverse  stresses  in  shocked  missile  material  (psi) 

x  y  r 
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oz  Axial  stress  in  missile  due  to  impact  shock  (psi) 

t  Time  parameter  which  describes  the  expansion  wave  profile  curve  on  the 

missile's  rigid  face  (sec.) 

<t>  Angle  between  deformed  missile  face  and  the  line  between  the  expansion 

wave  and  water  surface  contact  points  in  the  plane  of  symmetry  (rad.) 

^  Angle  between  wetting  side  of  missile  and  the  line  between  the 

expansion  wave  and  water  surface  contact  points  in  the  plane  of 
symmetry  (rad.) 

u  Deformation  angle  of  missile  face  from  its  rigid  profile  (rad.) 
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APPENDIX  A 

'IMPACT'  PROGRAM  LISTING 


82/09/02.  14.45.22. 

PRBGRAM  IMPACT 

00100  PR8GRAM  IMPACT  (INPUT, IUTPUT) 

00110  COMMON  /TABLE/  N,D 

00120  CONMON  /PHYS/  RHO ,1/1 , ALPHA, RH0K, AM, EZ, OMEGA, U, A, BETA, G,P 
00130  REAL  NU 
001  ADC 

001 50C  CONSTANT  PHYSICAL  DATA  ENTRY 

00160C 

001 70C  RICE  *  WALSH  SHOCK  HAVE  SPEED  DATA 
00180  DIMENSION  D(2,30) 

00190  DATA  N/30/ 

00200  DATA  (D<1, I), 1*1, 30)/-72.3E3,G. ,?I.5E3 , 1 45.0E3 ,21 7.5E3 ,290.0E3,342.5E3, 
0021 0+435. 0E3, 580. 0E3, 725. 0  E3 ,  Ej  70 .  DE3  f  1 . 0 1 5  E  6  P*  1 . 1 AOEA,  1  . 305E4, 1  .450E4, 

002 20+1. 595 E4,1.740E4,1  . 885E4,2.030E4,2.t?5E4,2.445E4,2.755E4,2.90QE6, 
OO230+3.O45E6,3.335E4,3.425E4,4.35OE4,5.0?5E4,5.B0OE4,4.525E4/ 

00240  DATA  ( D ( 2, 1 ) , I *1 , 30 1/25 . 98E3 ,58. 39E3, 79 .03E3, 92 . 40E3 , 1 03 . 4E3 , 1 1 3 .0E3 , 
00250+121 .4E3,128.9E3,142.2E3,153.9E3,144.3E3,1?3.BE3,182.4E3,190.9E3, 

00240  +  198. 4E3, 205. 8E3, 212. BE3, 219. 3E 3,225. 6E3, 231. 6E 3, 242. 5E3, 2  (53. 1E3, 

00270+258. 5E3, 243. 0E3, 273. 3E3, 281. 5E3, 301. 9E3, 320. 1E3, 334. 7E3, 352. 1E3/ 

00280  DATA  PI/3 . 141 6/,RH«/93.41 E-4/ ,C1 /58.39E3/ 

00290  PRINT  1 

00300  1  FORMAT  (///IBLIQUE  »r<  NIRMAL  WATER  ENTRY  IMPACT  OF  FLAT  FACED  CYLINDRIC 
00310+AL  MIS8LE6  (CONST.  ENTRY  VEL.)') 

00320C 

00330C  ENTRY  OF  INDEPENDENT  VARIABLES 

00340C 

00350  PRINT  2 

00340  2  FORMAT  (//ENTER  RADIUS  IF  MISSILE,  IN  «') 

00370  READ  *,  R 
00380  PRINT  3 

00390  3  FORMAT  (/.'IMPACT  VELOCITY  OF  MISSILE,  FT/SEC.  =') 

00400  REAP  +,  92 
00410  91*12. +92 
(10420  PRINT  4 

00430  4  FIRMAT  </, 'ANGLE  OF  ENTRY  FRIM  VERTICAL,  PEG.  *') 

00440  REAP  *,  ALPHA  1 
00450  ALPHA«ALPHA1*PI/I80. 

00440  PRINT  7 

00470  7  FIRMAT  (/, 'ENTER  HEIGHT  DENSITY  IF  MISSILE,  LBUN+  +  3  •  ■' ) 
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00480  READ  +  ,  RHIM1 
00490  RH0M-RH0M1 /38A.07 
00500  PRINT  8 

00510  B  FORMAT  ( /, '  YOUNGS  ELASTIC  MODULUS,  PSI  «') 

00520  READ  *,  E 
00530  PRINT  9 

00540  9  FORMAT  (/, 'POISSONS  RATIO  FOR  MISSILE  •') 

00550  READ  *,  NU 
00540  PRINT  11 

00570  11  FORMAT  </,  -'ENTER  THE  NUMBER  OF  TIME  INTERVALS  DESIRED  =") 

00380  READ  *,  M 
00590C 

Q0400C  OUTPUT  TABLE  HEADER 

00410C 

00420  PRINT  12 

00430  12  FORMAT  (//,17X, 'WATER  ENTRY  IMPACT  OF  FLAT  FACED  CYLINDRICAL  MISSILES  ( 
00440+C0N8T .  ENTRY  VEL.  >',/,43X, 'SHOCK  DYNAMICS  RESPONSE") 

00450  PRINT  13,  R,RH0M1 , ALPHA1 ,E, V2,NU 

00440  13  FORMAT  </ , 7X, 'MISSILE  RADIUS  *  ',F4.3,'  IN' , 13X, 'WEIGHT  DENSITY  * 
00470+F5.3,'  LB/IN**3',/,7X, 'ENTRY  ANGLE  *  ',F4.3('  DEG 1 4X, 'YOUNGS  MODULUS  ■ 
00480+', E9. 3,'  PSI', /,7X, 'ENTRY  VELOCITY  -  ',F5.0,'  FT/SEC. ',9X, 'POISSONS  RATIO 
00490+a  ' ,F5.3) 

00700C 

0071 OC  COMPUTATION  OF  CONSTANT  SHOCK  DYNAMICS  RESULTS 

00720C 

00730C  ELASTIC  WAVE  SPEED  IN  MISSILE 

00740  AM-SQRT(E*(1.-NU)/(RH0M*(1.+NU)+(1 .-2.*NU>>> 

00750C 

00740C  SHOCK  PRESSURE  ITERATIVE  COMPUTATION 

00770C  INITIAL  TRIAL  SHOCK  PRESSURE,  PI 

00780  J«1 

00790  P1»RH0*CI*V1 
00800  CALL  PRESS < P 1 ,01 ) 

00810  P2«P1+0.I*P1 
00820  15  CALL  PRE88(P2,G2) 

00830  CALL  SECANT  (0.0001 ,30, J, PI ,GI ,P2,G2) 

00840  IF  (J.GT.G)  G0  TC  15 
00850C 

OOB40C  CONVERGENT  VALUE  IF  SHICK  PRESSURE  Tf  +,-  0.01  PER  CENT,  P 
00870  CALL  PRE8S(P2,G) 

00880C 

00890C  EXPANSIIN  WAVE  SPEED  ACRISS  FACE  IF  MISSILE,  C 

00900  DELTA-PI/2. +0MEOA 

00910  PHI* AO  INI (V1-EZ*AM)+SIN(DELTA)/C1 ) 

00920  P8I«PI-(DELTA+PHI> 

00930  C»C1*8IN(PS1 ) /SIN (DELTA )*CIS(0NEGA) 

00940C 

00950C  TIME  AT  FULL  WETTING  ENTRY 
00940  TE*2.+R*TAN(ALPHA)/V1 
00970C 

00980C  TIME  DURATION  IF  IMPACT  SHICK 
00990  TI»R*TAN(ALPHA)/V1+R/C 
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oioooc 

OIOIOC  PRINT  #UT  AF  CONSTANT  SHACK  DYNAMICS  RESULTS 
01020  AMEGA 1 *AhEGA* 1 80 . /PI 
01030  BETA) »BETA*1B0./PI 

01040  PRINT  5,  AMfE2,0MEOA1fU,AfBETA1,Q(P,CfTE,TI 

01050  5  FiRMAT  (//SHACK  WAVE  SPEED  IN  MISSILE  *',E10.4/  IN/SEC/,/, 

01060+ 'ELASTIC  STRAIN  AN  MISSILE  FACE  »',£10.4/  IN/IN',/, 

01070+'DEFARMATIIN  ANGLE  AF  MISSILE  FACE  *',E10.4/  DEQ 
01080+  'WATER  PARTICLE  VELACITY  IN  SHACK  «',E10.4/  IN/SEC.  V, 

01090+'SHACK  UAVE  SPEED  IN  WATER  «',E10.4/  IN/SEC/,/. 

Oil OO+'MACH  ANGLE  AF  SHACK  UAVE  IN  WATER  =  ',E10.4/  DEG. ',/, 

0111 0+ "WATER  PARTICLE  SPEED  NARMAL  TA  SHACK  UAVE  »',Et0.4/  IN/SEC. 

01 120+'RANKINE-HUGANIAT  SHACK  PRESSURE  »',E10.4/  P SI',/, 

01130+'EXPANSIAN  WAVE  SPEED  ACRA9S  MISSILE  FACE  »',E10.4/  IN/SEC.',/, 
OIHO+'TIHE  AT  FULL  WETTING  ENTRY  «',E10.4/  SEC.'/, 

01150+'TIHE  DURATIAN  AF  IMPACT  SHACK  «',E10.4/  SEC/) 

01160C 

01 1 70C  CAMPUTATIAN  AF  TIME  DEPENDENT  RESULTS 

01180C 

01190  PRINT  A 

01200  A  FARHAT  </, AX,  'TIME' , 10X /SHACK  AREA' ,  4X/SHACK  FARCE', 3X, 
01210+'IHPULSE',/,6X/SEC/,10X/IN**2',9X/LB',l2X/LB-SEC/,/> 

O1220C 

01230C  TIME  STEP  INTERVAL,  DELTA  T 

01240  H«TI/<2.*N) 

01250C 

01260C  INITIALIZE  IMPULSE  ' 

01270  FT*0. 

01290  L“2*M+1 
01290  DA  1A  1*1 ,L 
01300  T=(I-1  )*H 
01 31 OC 

01320C  P^ACK  PRESSURE  AREA  DETERMINATION,  S 

01330  IF  (ALPHA. GT.O.)  G0  TA  10 

01340  S*PI*(R-C+T)+*2 

01350  OA  TA  19 

013AOC 

01370C  WETTING  LINE  SPEED 

01380  18  V=V1/TAN(ALPHA) 

01390C 

01  40 0 C  (JETTING  LINE  SPEED  CHECK 
OHIO  IF  (V.GT.C)  GO  TA  10 
0H20  PRINT  20,  V,C 

01430  20  FARMAT  (/.//THE  (JETTING  LINE  SPEED' ,E1 0. 4  /  IN/SEC.',  'IS  LESS  THAN  THE 
01440+EXPANBIAN  WAVE  SPEED',E10.4/  IN/SEC. ", //NA  IMPACT  SHACK  ACCURS') 

DH50  STAP 

01460  10  CALL  AREA  (R,V,C,T,200,S1 ) 

01470  S=2.+S1 
CM480C 

01490C  IMPACT  SHACK  FARCE 

01500  19  F  =  P*S 
01 1 5 1 0C 

01520C  INTEGRATE  FARCE  T«  IHPULSE 

01530  CALL  SI  MR  (I,H,F/T) 
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01 540C 

01550C  PRINT  0UT  IF  SHICK  DYNAMICS  RESULTS 

01540  IF  < I .EQ.2*( 1/2) )  60  T0  14 

01570  PRINT  17,  T,S,F,FT 

01500  14  CONTINUE 

01590  17  FIRNAT  (4<4X rE 1 0. 4 > ) 

01400  END 
01 41 OC 

01 420C  RANKINE-HUG0NIIT  SHICK  PRESSURE  COMPUTATION 

O1430C 

01440  SUBROUTINE  PRESS  <P ,Q > 

01450  C0HM0N  /PHYS/  RH0,V1 , AL PHA , RH0H , AM , E Z , BHEOA, U , A, BETA , Qf PS 
014400 

01470C  AXIAL  STRAIN  IN  HISSILE  UNDER  SH0CK  PRESSURE 
01480  EZ*P/(RHIM*AH*AH) 

01490C 

0170DC  DEFIRMATIIN  ANGLE  0F  NISSILE  FACE  FRIM  RIGID  FACE 
01710  IHEGA*ATAN< EZ*AM* TAN ( ALPHA) /VI > 

01720C 

01730C  UATER  PARTICLE  VELICITY 

01740  U=(V1-EZ*AM>*CIS<0HEGA) 

01 7500 

O1740C  UATER  SH0CK  SPEED  FRIM  TRIAL  PRESSURE  (RICE  3  UALSH) 

01770  CALL  ITABLE  < P, A> 

01780C 

Q1790C  SURFACE  DISTURBANCE  SPEED  CHECK 

01B00  IF  (VI /A. 0T. SIN< ALPHA ) )  80  Tf  2 
01810  D8=V1 /S I N < ALPHA ) 

01820  PRINT  1,  D5, A 

01830  1  F0RNAT  (///SURFACE  DISTURBANCE  SPEED , E 1 0 . 4 , '  IN/SEC.  IS  LESS  THAN  THE 
O1B4O+SH0CK  WAVE  SPEED" ,E10.4, "  IN/SEC. ",/, "N0  IMPACT  SH0CK  0CCURS . " > 

01850  ST0P 
01840C 

01870C  UATER  SH0CK  MACH  ANGLE 

01880  2  BETA«ASIN(A*8IN(ALPHA)/V1) 

01890C 

01900C  C0NP0NENT  IF  UATER  PARTICLE  VELICITY  N0RMAL  T0  3H0CK  FRINT 
01910  G“U+C0S< BETA—  <ALPHA— BMEGA) > 

01920C 

01930C  RANKINE-HUOINI0T  SHICK  PRESSURE i  0BLIGUE  SHICKS 

01940  P3=RH0*A*Q 

01950C 

01960C  SHICK  PRESSURE  ERRIR  FUNCTI0N 

01970  G-P-PS 
01980  RETURN 
01990  END 

02000C  -  *• 

0201  OC  4  PUNT  LAGRAN8IAN  INTERP0LATI0N  0F  THE 

Q2020C  RICE  1  UALSH  SHICK  PRESSURE  3  UAVE  SPEED  DATA  TABLE 

02030C 

02040  SUBR0UTINE  ITABLE  (X,Y) 

02050  C0HN0N  /TABLE/  N,T 
02040  DIHENSIIN  T(2,1) 
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02070  H*N-2 
02080  P0  1  1=2, H 

02090  IF  (X.GT.T(1,I).ANB.X.LE.T(t,I*1)>  G 0  TO  3 
02100  1  CONTINUE 
02110  PRINT  2,  X 

02120  2  FORMAT  </,'X  * " , E 1  □ . 4 ,  '  LIES  OUTSIPE  THE  INNER  RANGE  IF  THE  TABLE1' ) 

02130  STOP 

02140  3  XI =T< 1 , 1-1 ) 

02150  V1=T(2,I-1) 

02160  X2=T(1,I) 

02170  Y2=T (2,1) 

02180  X3*T (1,1+1) 

02190  Y3“T (2,1  +  1 ) 

02200  X4=T ( 1 ,1+2) 

02210  Y4«T(2, 1+2 ) 

02220  M>(X1-X2)*(Xt*X3)*(X1-X4) 

02230  D2* < X2-X 1 > * (X2-X3 )*<X2-X4> 

02240  P3*<  X3-X1 )*(X3-X2)*(X3-X4 ) 

02250  D4*(X4-X1)*<X4-X2)*(X4-X3) 

02260  V«(X-X2)*(X-X3)*(X-XO*V1/»1 
02270  Y«Y4(X-X1 >*(X-X3)*(X-X4)*V2/B2 
02280  Y»Y+(X-X1 )*<X-X2)*(X-X4)*Y3/B3 
02290  Y«r+(X-X1 )*<X-X2)*(X-X3>*Y4/D4 
02300  RETURN 
02310  ENP 
02320C 

023300  9ECANT  ITERATION  OF  SHOCK  PRESSURE  Tl  YIELD  A  RIOT 

023400  IF  THE  PRESSURE  ERROR  FUNCTION  G<P)«0. 

02350C 

U2360  SUBROUTINE  SECANT(E,N,I,X1 ,F1 ,X2,F2) 

02370  IF  (ABS(F2-F1).EQ.O.)  00  TO  1 
02380  X3»X2-(X2-X1)*F2/<F2-F1> 

02390  X1-X2 
02400  F 1=F2 
02410  X2»X3 
02420  1*1+1 

02430  1  IF  (APS((X2-X1)/Xi).LE.c)  G0  TO  2 

02440  IF  (I.GT.N)  GO  Tl  3 

02450  RETURN 

02440  2  I*-1 

02470  RETURN 

02480  3  PRINT  4,  N,X1 ,F1 ,X2,F2 

02490  4  FORMAT  U, -SECANT  METHOD  POES  NOT  C0NVERBE  FOR  N  “'',13,/, 

02500+  'XI  = ',E12.6,4X, 'FI  ='fE12.6,/'X2  *',E12.4,4X,  'F2  "=',E12.6) 

02510  STOP 
02520  ENP 
02530C 

O2540C  HALF  AREA  BETWEEN  WETTING  LINE  «  EXPANSION  WAUE  FRONT  PROFILE 

02550C 

02560  SUBROUTINE  AREA  <R, V,C,T,N,A> 

02570  A  =  0 . 

02580  IF  (T.LE.0..0R.T.0E.R/V+R/C-1E-9)  GO  Tl  4 
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02590  X1»0. 

Q26Q0  X2*Q. 

02410  Y1«0. 

02420  Y2-0. 

02430  T1»T 

02440  IF  iT.QT.2.*R/V)  T1»2.»R/V 
02450  M 3 M+ 1 
02440  00  1  1  =  1, H 
0  2  4  ?  0  C 

O2480C  TIME  PARAMETER  F0R  PROFILE  CURVE,  0<TAU<T1 
02490  TAU- <1-1  >MT1/N) 

02700  XI =R-V*TAU 
02710  ETA«80RT<R*R-XI*XI> 

02720  RH0=C*<T-TAU) 

02730  DELTA“IC/V)+RHI 

02740  RT-SQRT(R*R*RHI*RHf-ETA*ETA*DELTA«DELTA> 

Q2750C 

02740C  COORDINATES  0F  PROFILE  CURVE*  X(TAU) , Y(TAU) 

02770  X»XMETA+ETA+DELTA+X1*RT)/(R*R> 

02780  Y=ETA*<X+DELTA)/XI 
02790C 

028000  SIMP90N8  RULE  INTEGRATION  FOR  UNEQUAL  INTERVALS  ON  X 

02810  IF  (I.6T.1)  80  TO  2 

02820  X1-X 

02830  Y 1  * Y 

02840  GO  TO  1 

02850  2  IF  ( I  ,GT.2+(I/2) )  GO  TO  3 
02840  X2*X 

02870  Y2-Y  , 

02880  GO  TO  1 

02890  3  P1=(X1**3-X*+3>/3. 

02900  P2=(X1*X1“X*X)/2. 

02910  P3-X1 -X 

02920  A»A+ ( PI  - ( X2+^)+P2+X2*X*P3  )*Y1/(<X1-X2)*(X1-X>) 

02930  A«A+tP1-(X1+X)*P2+X1*X*P3)*Y2/((X2-X1 >*(X2-X>> 

02940  A=A+ ( PI- (X1+X2) *P2+XI *X2*P3)*Y/( { X-X 1 >*(X-X2>  > 

02950  XI -X 
02960  Y 1  * Y 
02970  1  CONTINUE 
02980  4  RETURN 
02990  END 
03000C 

0301 OC  SIMPSON'S  1/3  RULE  INTEGRATION 

03020C 

03030  SUBROUTINE  SIMP  <I,H,Y,A> 

03040  IF  (I.QT.1)  GO  TO  1 
03050  Y1 -Y 
03060  RETURN 

03070  I  IF  (I. OT. 2*11/2))  GO  TO  2 
03080  Y2«Y 
03090  RETURN 

03100  2  A»A+(Y1*4.*Y2+Y)*H/3. 

03110  Y 1 »Y 
03120  RETURN 
03130  END 
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APPENDIX  B 

LISTING  OF  EXAMPLE  OUTPUT 


run 

82/09/02.  16.52.03. 
PROGRAM  IMPACT 


OBLIQUE  OR  NORMAL  HATER  ENTRY  IMPACT  OF  FLAT  FACED  CYLINDRICAL  MISSLES  (CONST.  ENTRY  VEL.) 

ENTER  RADIUS  OF  MISSILE,  IN  * 

T  6 

IMPACT  VELOCITY  OF  MISSILE,  FT/SEC.  » 

?  1-500 

ANOLE  IF  ENTRY  FROM  VERTICAL,  DEO.  » 

*  5 

ENTER  UEIGHT  DENSITY  OF  MISSILE,  LB/IN*+3  * 

T  0.1 

YOUNSS  ELASTIC  MODULUS,  F'SI  * 

*  10*0  ■ 

POISSONS  RATIO  FOR  MISSILE  * 

?  0.333 

ENTER  THE  NUMBER  OF  TIME  INTERVALS  DESIRED  » 

?  50 


HATER  ENTRY  IMPACT  OF  FLAT  FACED  CYLINDRICAL  MISSILES  (CONST.  ENTRY  VEL.) 

SHOCK  DYNAMICS  RESPONSE 


MISSILE  RADIUS  ‘  6.000  IN 

ENTRY  ANGLE  *  5.000  DEG. 

ENTRY  VELOCITY  =  1500.  FT/SEC. 


HEIGHT  DENSITY 
YOUNGS  NODULUS 
POISSONS  RATIO 


.100  LB/IN*#3 
.IOOE+OB  PSI 
.333 


SHOCK  WAVE  SPEED  IN  MISSILE  «  .2405E+06  IN/SEC. 
ELASTIC  STRAIN  ON  MISSILE  FACE  *  .S301E-02  IN/IN 
DEFORMATION  ANGLE  OF  MISSILE  FACE  *  .5560E+00  DEO. 

' UATER  PARTICLE  VELOCITY  IN  SHOCK  =  .1600E+05  IN/SEC. 
SHOCK  WAVE  SPEED  IN  HATER  =  .8920E+05  IN/SEC. 

NACH  ANGLE  OF  SHOCK  HAVE  IN  WATER  »  .2559E+02  DES . 
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UATER  PARTICLE  SPEED  N0RMAL  T0  SH8CK  UAUE  *  .1493E+05  IN/SEC. 
RANK INE-HUQf NI0T  SH0CK  PRESSURE  ■«  .1244E+04  PSI 
EXPANSION  WAVE  SPEED  ACR0SS  MISSILE  FACE  =  .340DE+05  IN/SEC. 
TIME  AT  FULL  UETTIN0  ENTRY  =>  .SB33E-04  SEC. 

TIME  DURATION  0F  IMPACT  SHOCK  -  .13A3E-03  SEC. 


TIME 

SEC. 

0. 

.  2724E-05 

.  5453E-05 
.8I79E-05 
.  1 091 E“04 
.  1 3A3E-04 
.  1 434E-04 
.  1 908E-O4 

j  .2181E-04 

1  .2454E-04 

.  272AE-04 

!  .2999E-04 

i  .3272E-04 

.3544E-04 

1  .3817E-04 

.4089E-04 

>  .  43A2E-04 

|  , 4435E-04 

'  .4907E-04 

J  .51 80E-04 

.3433E-04 

;  .5723E-04 

S  .599BE-04 

;  .A270E-04 

•4543E-04 

i  .681 6E-04 

!  .7088E-04 

.73A1E-04 

'  .7A34E-04 

|  .7904E-04 

i  .81 79E-04 

.B431E-04 

•8724E-04 

.B997E-04 

.9249E-04 

1  .9542E-04 

.9815E-04 
.1009E-03 
. 1036E-03 

1  . 1043E-03 

. 1 091 E— 03 
.111 8E-03 
. 1 1 45E-03 


SHICK  AREA 

SHOCK  FORCE 

IMPULSE 

I M  *  2 

LB 

LB-SEC. 

■ 

0. 

0. 

.  1 387E+0 1 

. 1 725E+04 

.  1 902E+00 

. 3B50E+01 

.4788E+04 

. 1 Q37E+01 

. A935E+0 1 

. 8425E+04 

. 2870E+O1 

. 1 Q44E+Q2 

.1 301E+07 

.5809E+01 

.  1 43 1 E+Q  2 

. 1 780E+07 

. 1 OOOE+02 

.  1 84  0E+0  2 

. 2  289E+07 

.1354E+02 

, 2245E+02 

.281 7E+07 

, 2250E+02 

.27Q1E+Q2 

.3359E407 

.  3092E+02 

, 31 40E+02 

. 3905E+0? 

.4082E+02 

. 3578E+02 

•4430E+07 

. 5221 E+02 

.401 0E+02 

. 4987E+07 

.  6507E+02 

.4429E+02 

.3509E+07 

.7938E+02 

. 4831E+02 

.40O8E+O7 

.  9509E+02 

.5209E+02 

.4478E+07 

.1121E+Q3 

.5537E+02 

. 4910E+07 

.  1 304E+O3 

. 5844E+02 

.729AE+0? 

.  1 498E+03 

. 41 30E+02 

. 7A23E+07 

.1701 E+03 

. A335E+02 

. 7879E+0? 

.  1 91 3E+03 

. 444BE+02 

.8044E+07 

.2130E+03 

.4504E402 

•8092E+0? 

.  2350E+03 

. 6405E+02 

.7945E+07 

. 2549E+03 

.4043E+02  ' 

.7514E+07 

.2781 E+03 

.543 1 E+02 

.7002E+07 

.2979E+03 

.  5232E+D2 

.4507E+07 

.  31 43E+03 

.4B47E+02 

.4028E+C7 

•3334E+03 

.4478E+02 

•5549E+07 

.  3492E+03 

.41 22E+02 

.5I24E+07 

.  3438E+03 

.3782E+02 

. 4704E+07 

.  3772E+03 

.  3457E+02 

.4299E+07 

.  389 4 E+03 

.31 44E+02 

.3913E+07 

.  4004E+03 

.  2830E+02 

.3543E+07 

•4108E+03 

.  2548E  +  02 

.3193E+07 

.  4200E+03 

.  2302E+02 

. 2843E+07 

.  4282E+03 

.  2048E+02 

•2544E+07 

.4334E+03 

.  181 3E+02 

•2255E+07 

.  4421 E+03 

.  1 588E+02 

. 1 975E+D7 

.4479E+03 

.  I380E+02 

•1714E+07 

.  4529E+03 

.  1 1 84E  +  02 

. 1 473E+Q7 

.4573E+03 

.  1003E+02 

.1250E+0? 

.  441 0E+03 

.8417E+01 

. 1047E+07 

.  4641 E+03 

.491 4E+01 

.8401 E+04 

. 4467E+03 

.3551 E+01 

.49048+06 

.  448QE+03 
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. 1 172E-03 

•4343E+01 

.5401 E+OA 

.4705E+O3 

. 1 200E-03 

.3288E+01 

. 4088E+O4 

.  471 8E+03 

.  1 227E-03 

.2344E+01 

.2941E+04 

.4727E+03 

.  1  254  E -03 

. 1 594E+01 

. 1 982E+04 

.  4734E+03 

. 1 281 E-D3 

.9492E+00 

. 1 205E+04 

.4738E+03 

. 1 309E-03 

.4895E+0Q 

.A087E+O5 

.  474 1 E+03 

. 1 33AE-03 

. 1 572E+00 

. 1 9S5E+05 

.  4742E+03 

. 1 363E-03 

0. 

0. 

.4742E+03 

8RU  10.185  UNTS. 

RUN  C0HPLETE. 
by* 

249  LIO  0FF  16.55.32. 

249  SRU  10.382  UMTS. 
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